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1. Introduction 



The experimental efforts at RHIC, Q liave provided a novel window in the physics of the 
strong interactions. The consensus on the existing data is that shortly after the collision, a 
ball of quark-gluon plasma (QGP) forms that is at thermal equilibrium, and subsequently 
expands until its temperature falls below the QCD transition (or crossover) where it finally 
hadronizes. Relativistic hydrodynamics describes very well the QGP with a shear- 
viscosity to entropy density ratio close to that of A/" = 4 SYM, The QGP is at 
strong coupling, and it necessitates a treatment beyond perturbative QCD approaches, Q. 
Moreover, although the shear viscosity from A/" = 4 seems to be close to that "measured" 
by experiment, lattice data indicate that in the relevant RHIC range 1 < f < 3 the QGP 
seems not to be a fully conformal fluid. Therefore the bulk viscosity may play a role near 
the phase transition ||5|, |^. The lattice techniques have been successfully used to study 
the thermal behavior of QCD, however they are not easily extended to the computation 
of hydrodynamic quantities. They can be used however, together with parametrizations 
of correlators in order to pin down parameters Q. On the other hand, approaches based 
on holography have the potential to address directly the real-time strong coupling physics 
relevant for experiment. 

In the bottom- up holographic model of AdS/QCD the bulk viscosity is zero as 
conformal invariance is essentially not broken (the stress tensor is traceless) . In the soft- wall 
model [^, no reliable calculation can be done for glue correlators and therefore transport 
coefficients are ill-defined. Similar remarks hold for other phenomenologically interesting 
observables as the drag force and the jet quenching parameter |^, ^ [T^ . 

Top-down holographic models of QCD displaying all relevant features of the theory 
have been difficult to obtain. Bottom- up models based on AdS slices [13] have given some 
insights mostly in the meson sector, |^ but necessarily lack many important holographic 
features of QCD. A hybrid approach has been advocated |l5| combining features of 
bottom-up and top-down models. An similar approach was proposed independently in [p!6| . 
Such an approach, called Improved Holographic QCD (or IHQCD for short) is essentially 
a five-dimensional dilaton-gravity system with a non-trivial dilaton potential. Flavor can 
be eventually added in the form of Nf space-time filling D4 — DA brane pairs, supporting 
U{Nf)L X U{Nf)R gauge fields and a bi-fundamental scalar [l?]. The UV asymptotics 
of the potential are fixed by QCD perturbation theory, while the IR asymptotics of the 
potential can be fixed by confinement and linear glueball asymptotics. An analysis of the 
finite temperature behavior |l^ has shown that the phase structure is exactly what 
one would expect from YM. A potential with a single free parameter tuned to match the 
zero temperature glueball spectrum was able to agree with the thermodynamic behavior 
of glue to a good degree, |18]. Similar results, but with somewhat different potentials were 
also obtained in [16, 2C] 



In |18|, |19|, [ig] it was shown that Einstein-dilaton gravity with a strictly monotonic 
dilaton potential that grows sufficiently fast, generically shares the same phase structure 
and thermodynamics of finite-temperature pure Yang-Mills theory at large Nc- There is a 
deconfinement phase transition (dual to a Hawking-Page phase transition between a black 
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hole and thermal gas background on the gravity side), which is generically first order. The 
latent heat scales as N"^. In the deconfined gluon-plasma phase, the free energy slowly 
approaches that of a free gluon gas at high temperature, and the speed of sound starts 
from a small value at Tc and approaches the conformal value = 1/3 as the temperature 
increases. The deviation from conformal invariance is strongest at Tc, and is signaled by the 
presence of a non-trivial gluon condensate, which on the gravity side emerges as a deviation 
of the scalar solution that behaves asymptotically as close to the UV boundary. In the 
CP-violating sector, the topological vacuum density TrFF has zero expectation value in 
the deconfined phase, in agreement with lattice results [^] and large- iVc expectations. 

The analysis performed in [^] was completely general and did not rely on any specific 
form of the dilaton potential V{\). A detailed analysis of an explicit model in ||2^ shows 
that the thermodynamics matches quantitatively the thermodynamics of pure Yang-Mills 
theory. The (dimensionless) free energy, entropy density, latent heat and speed of sound, 
obtained on the gravity side by numerical integration of the 5D field equations, can be 
compared with the corresponding quantities, calculated on the lattice for pure Yang-Mills 
at finite-T, resulting in excellent agreement, for the temperature range that is accessible 
by lattice techniques. The same model also shows a good agreement with the lattice 
calculation of glueball mass ratios at zero temperature, and the value of the deconfining 
critical temperature (in units of the lowest glueball mass) is also in good agreement with 
the lattice results. 

In short, the model we present gives a good phenomenological holographic description 
of most static properties^ (spectrum and equilibrium thermodynamics) of large- A'^c pure 
Yang-Mills, as computed on the lattice, for energies up to several times Tc. Thus it con- 
stitutes a good starting point for the computation of dynamical observables in a realistic 
holographic dual to QCD (as opposed to e.g. M = A SYM), such as transport coefficients 
and other hydrodynamic properties that are not easily accessible by lattice techniques, at 
energies and temperatures relevant for relativistic heavy-ion collision experiments. We will 
report on such a calculation in the near future. 

The vacuum solution in this model is described in terms of two basic bulk fields, the 
metric and the dilaton. These are not the only bulk fields however, as the bulk theory is 
expected to have an a priori infinite number of fields, dual to all possible YM operators. In 
particular we know from the string theory side that there are a few other low mass fields, 
namely the RR axion (dual to the QCD 0-angle) the NSNS and RR two forms B2 and C2 
as well as other higher-level fields. With the exception of the RR axion, such fields are dual 
to higher-dimension and/or higher-spin operators of YM. Again, with the exception of the 
RR axion, they are not expected to play an important role into the structure of the vacuum 
and this is why we neglect them when we solve the equations of motion. However, they 
are going to generate several new towers of glueball states beyond those that we discuss 
in this paper (namely the glueballs associated to dilaton fluctuations, 2++ glueballs 
associated to graviton fluctuations and ^ glueballs associated to RR axion fluctuations) . 
Such fields can be included in the effective action and the associated glueball spectra 

^ There are very few observables also that are not in agreement with YM. They are discussed in detail 
in ||. 
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calculated. Since we do not know the detailed structure of the associated string theory, 
their effective action will depend on more semi-phenomenological functions like Z(\) in 
( |5.9| ). These functions can again be determined in a way similar to Z{X). In particular 
including the B2 and C2 field will provide 1"' glueballs among others. Fields with spin 
greater than 2 are necessarily stringy in origin. We will not deal further with extra fields, 
like B2 and C2 and other as they are not particularly relevant for the purposes of this 
model, namely the study of finite temperature physics in the deconfined case. We will only 
consider the axion, as its physics is related to the CP-odd sector of YM with an obvious 
phenomenological importance. 

It is well documented that string theory duals of YM must have strong curvatures in 
the UV regime. This has been explained in detail in where it was also argued, that 
although the asymptotic AdS boundary geometry is due to the curvature non-linearities 
of the associated string theory, the inwards geometry is perturbative around AdS, with 
logarithmic corrections, generating the YM perturbation theory. The present model is 
constructed so that it takes the asymptotic AdS geometry for granted, by introducing the 
associated vacuum energy by hand, and simulates the perturbative YM expansion by an 
appropriate dilaton potential. In the IR, we do not expect strong curvatures in the string 
frame, and indeed the preferred backgrounds have this property. In this sense the model 
contains in itself the relevant expected effects that should arise from strong curvatures in 
all regimes. These issues have been explained in [|lj] and in more detail in [^]. 

A different and interesting direction is the use of such models to study the expansion 
of the plasma and the associated dynamics. Such a context is similar to what happens 
on cosmology, especially the one related to the Randall-Sundrum setup. Indeed in this 
case the expansion can be found by following the geodesic motion of probe branes in the 
relevant background, |23|-[^]. This generalizes to more complicated backgrounds, |26|, 
like the ones studied here. 

Once we have a holographic model we trust, we should calculate observables , like 
transport coefficients that are hard to calculate on the lattice. A first class of transport 
coefficients are viscosity coefficients.^ A general fluid is characterized by two viscosity 
coefficients, the shear r] and the bulk viscosity C. The shear viscosity in strongly coupled 
theories described by gravity duals was shown to be universal, Q. In particular, the ratio 
r]/s, with s the entropy density, is equal to This is correlated to the universality of low- 
energy scattering of gravitons from black holes. It is also known that deviations from this 
value can only be generated by higher curvature terms that contain the Riemann tensor 
(as opposed to the Ricci tensor of the scalar curvature). In QCD, as the theory is strongly 
coupled in the temperature range Tc < T < 3Tc, we would expect that rj/s ^ Recent 
lattice calculations, [^] agree with this expectations although potential systematic errors 
in lattice calculations of transport coefficients can be large. 

Conformal invariance forces the bulk viscosity to vanish. Therefore the = 4 SYM 
plasma, being a conformal fluid, has vanishing bulk viscosity. QCD on the other hand is not 



^ These are the leading transport coefficients in the derivative expansion. There are subleading coefficients 
that have been calculated recently for TV = 4 SYM, |^^. However, at the present level of accuracy, they 
cannot affect substantially the comparison to experimental data, |H]. 
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a conformal theory. The classical theory is however conformally invariant and asymptotic 
freedom implies that conformal invariance is a good approximation in the UV. This would 
suggest that the bulk viscosity to entropy ratio is negligible at large temperatures. However 
it is not expected to be so in the IR: as mentioned earlier lattice data indicate that in the 
relevant RHIC range 1 < £ < 3 the QGP seems not to be a fully conformal fluid. Therefore 
the bulk viscosity may play a role near the phase transition. 

So far there have been two approaches that have calculated the bulk viscosity in 
YM/QCD, 30, 31, 32| and have both indicated that the bulk viscosity rises near the 
phase transition as naive expectation would suggest. The first used the method of sum 



rules in conjunction with input from Lattice thermodynamics, 3C, 31|. It suggested a 
dramatic rise of the bulk viscosity near Tc although the absolute normalization of the result 
is uncertain. The reason is that this method relies on an ansatz for the density associated 
with stress-tensor two point functions that are otherwise unknown. 

The second method |32| relies on a direct computation of the density at low frequency 
of the appropriate stress-tensor two-point function. As this computation is necessarily 
Euclidean, an analytic continuation is necessary. The values at a finite number of dis- 
crete Matsubara frequencies are not enough to analytically continue. An ansatz for the 
continuous density is also used here, which presents again a potentially large systematic 
uncertainty. 

Calculations in IHQCD support a rise of the bulk viscosity near Tc, but the values 
are much smaller than previously expected. Studies of how this affects hydrodynamics at 
RHIC, suggest that this implies a fall in radial and elliptic flow. 

Another class of interesting experimental observables is associated with quarks, and 
comes under the label of "jet quenching". Central to this is the expectation that an 
energetic quark will loose energy very fast in the quark-gluon plasma because of strong 
coupling. This has as a side effect that back-to back jets are suppressed. Moreover if a 
pair of energetic quarks is generated near the plasma boundary then one will exit fast the 
plasma and register as an energetic jet, while the other will thermalize and its identity will 
disappear. This has been clearly observed at RHIC and used to study the energy loss of 
quarks in the quark-gluon plasma. 

Heavy quarks are of extra importance, as their mass masks some low-energy strong 
interaction effects, and can be therefore cleaner probes of plasma energy loss. There are 
important electron observables at RHIC, |34] that can probe heavy-quark energy loss in 
the strongly coupled quark-gluon plasma. Such observables are also expected to play an 
important role in LHC ||35[| . 

A perturbative QCD approach to calculate the energy loss of a heavy quark in the 
plasma has been pursued by calculating radiative energy loss, \2>t\. However its application 
to the RHIC plasma has recently raised problems, based on comparison with data. A 
phenomenological coefficient used in such cases is known as the jet quenching coefficient q, 
and is defined as the rate of change of the average value of transverse momentum square 
of a probe. Current fits, [34, 37|, indicate that a value of order 10 GeV"^ / fm or more is 
needed to describe the data while perturbative approaches are trustworthy at much lower 
values. 
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Several attempts were made to compute quark energy loss in the holographic context, 
relevant for A/" = 4 SYM^. In some of them |39, ^] the jet-quenching coefficient q was 
calculated via its relationship to a light-like Wilson loop. Holography was then used to 
calculate the appropriate Wilson loop. The q obtained scales as a/A and as the third power 
of the temperature, 

9conformal = [f | V2X TT^T^ (1.1) 

A different approach chooses to compute the drag force acting a string whose UV end- 
point (representing an infinitely heavy quark) is forced to move with constant velocity v, 
|9|, H, im , in the context of A/" = 4 SYM plasma. The result for the drag force is 

-^conformal = TT"^ , =^ (1-2) 

^ vi — 

and is calculated by first studying the equilibrium configuration of the appropriate string 
world-sheet string and then calculating the momentum flowing down the string. This can 
be the starting point of a Langevin evolution system, as the process of energy loss has a 
stochastic character, as was first pointed out in |41] and more recently pursued in |42]-|48|. 

Such a system involves a classical force, that in this case is the drag force, and a stochas- 
tic noise that is taken to be Gaussian and which is characterized by a diffusion coefficient. 
There are two ingredients here that are novel. The first is that the Langevin evolution must 
be relativistic, as the quarks can be very energetic. Such relativistic systems have been 
described in the mathematical physics literature, |49] and have been used in phenomeno- 
logical analyses of heavy-ion data, p7[ |. They are known however to have peculiar behavior, 
since demanding an equilibrium relativistic Boltzmann distribution, provides an Einstein 
relation that is pathological at large temperatures. Second, the transverse and longitudinal 
diffusion coefficients are not the same, A first derivation of such Langevin dynamics 
from holography was given in This has been extended in [^ ] where the thermal- like 
noise was associated and interpreted in terms of the world-sheet horizon that develops on 
the probe string. 

Most of the transport properties mentioned above have been successfully computed in 
A/" = 4 SYM and a lot of debate is still waged as to how they can be applied to QCD in 
the appropriate temperature range, A holographic description of QCD has 

been elusive, and the best we have so far have been simple bottom up models. 

In the simplest bottom-up holographic model known as AdS/QCD [^, the bulk viscos- 
ity is zero as conformal invariance is essentially not broken (the stress tensor is traceless), 
and the drag force and jet quenching essentially retain their conformal values. 

In the soft- wall model [P, no reliable calculation can be done for glue correlators 
and therefore transport coefficients are ill-defined, as bulk equations of motion are not 
respected. Similar remarks hold for other phenomenologically interesting observables as 
the drag force and the jet quenching parameter. 

The shear viscosity of IHQCD is the same as that of = 4 SYM, as the model is 
a two derivative model. Although this is not a good approximation in the UV of QCD, 



Most are reviewed in 
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it is expected to be a good approximation in the energy range Tc < T < 5Tc. The bulk 
viscosity in IHQCD rises near the phase transition but ultimately stays slightly below the 
shear viscosity. There is a general holographic argument that any (large-N) gauge theory 
that confines color at zero temperature should have an increase in the bulk viscosity-to- 
entropy density ratio close to Tc. 

The drag force on heavy quarks, and the associated diffusion times, can be calculated 
and found to be momentum depended as anticipated from asymptotic freedom. Numerical 
values of diffusion times are in the region dictated by phenomenological analysis of heavy- 
ion data. The medium-induced corrections to the quark mass (needed for the diffusion time 
calculation) can be calculated, and they result in a mildly decreasing effective quark mass as 
a function of temperature. This is consistent with lattice results. Finally, the jet-quenching 
parameter can be calculated and found to be comparable at Tc to the one obtained by 
extrapolation from Af = 4 SYM. Its temperature dependence is however different and 
again refiects the effects of asymptotic freedom. 



2. The 5D model 



The holographic dual of large Nc Yang Mills theory, proposed in |14], is based on a five- 
dimensional Einstein-dilaton model, with the action^: 



^5 = -M^N^ 



R- -{d^f + V{^) 



(2.1) 



dM 



Here, Mp is the five-dimensional Planck scale and Nc is the number of colors. The last 
term is the Gibbons-Hawking term, with K being the extrinsic curvature of the boundary. 
The effective five-dimensional Newton constant is G5 = 1/ {IQirMpNc), and it is small in 
the large- A^c limit. 

Of the 5D coordinates {xj, r}j=o...3i Xi are identified with the 4D space-time coordi- 
nates, whereas the radial coordinate r roughly corresponds to the 4D RG scale. We identify 
A = with the running 't Hooft coupling At = NcQym^ priori unknown multi- 

plicative factor^, A = nXt. 

The dynamics is encoded in the dilaton potential^, V^(A). The small-A and large- 
A asymptotics of V{\) determine the solution in the UV and the IR of the geometry 
respectively. For a detailed but concise description of the UV and IR properties of the 
solutions the reader is referred to Section 2 of |ll9[| . Here we will only mention the most 
relevant information: 



^Similar models of Einstein-dilaton gravity were proposed independently in to describe the finite 
temperature physics of large Nc YM. They differ in the UV as the dilaton corresponds to a relevant operator 
instead of the marginal case we study here. The gauge coupling e* also asymptotes to a constant instead 
of zero in such models. 

^This relation is well motivated in the UV, although it may be modified at strong coupling (see Section 
The quantities we will calculate do not depend on the explicit relation between A and At. 

®With a slight abuse of notation we will denote V(A) the function V{^) expressed as a function of 

\ — * 
A = e . 
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1. For small A, V{X) is required to have a power-law expansion of the form: 



y(A) ~ ^{l + voX + viX^ + ...) 



A ^ . 



(2.2) 



The value at A = is constrained to be finite and positive, and sets the UV AdS scale 
£. The coefficients of the other terms in the expansion fix the /3-function coefficients 
for the running coupling X{E). If we identify the energy scale with the metric scale 
factor in the Einstein frame, as in p4|, we obtain: 



/3(A) = 



dX 



dlogE 



-boX^ - hX^ + ... 



9 



207 



4 256 " 

2. For large A, confinement and the absence of bad singularities^ require: 



(2.3) 



y(A) ~ A2'3(logA)^ A^oo, 



2/3 <Q < 2^/2/3, P arbitrary 
Q = 2/3, P>0 



(2.4) 



In particular, the values Q = 2/3, P = 1/2 reproduce an asymptotically- linear glue- 
ball spectrum, mH ~ 
in what follows. 



^ ~ n, besides confinement. We will restrict ourselves to this case 



In the large Nc limit, the canonical ensemble partition function of the model just 
described, can be approximated by a sum over saddle points, each given by a classical 
solution of the Einstein-dilaton field equations: 



^(/3):^e-'5^(« + e-'5^('') + ... 



(2.5) 



where Si are the euclidean actions evaluated on each classical solution with a fixed temper- 
ature T = 1//3, i.e. with euclidean time compactified on a circle of length (3. There are two 
possible types of Euclidean solutions which preserve 3-dimensional rotational invariance. 
In conformal coordinates these are: 



1. Thermal gas solution, 

ds^ = bl{r) [dr^ + dt^ + dxmdx"^) , $ = $„(r), 
with r S (0, oo) for the values of P and Q we are using; 

2. Black-hole solutions, 



(2.6) 



ds' = b{ry 



fir. 



+ f{r)dr + dxjndx' 



$ = <I>(r), 



(2.7) 



with r G (0, r/i), such that /(O) = 1, and f{rh) = 0. 



^We call "bad singularities" those that do not have a well defined spectral problem for the fluctuations 
without imposing extra boundary conditions. 
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In both cases Euclidean time is periodic with period /3o and /3 respectively for the thermal 
gas and black-hole solution, and 3-space is taken to be a torus with volume V30 and V3 
respectively, so that the black-hole mass and entropy are finite^. 

The black holes are dual to a deconfined phase, since the string tension vanishes at 
the horizon, and the Polyakov loop has non- vanishing expectation value (|5^, On the 
other hand, the thermal gas background is confining. 

The thermodynamics of the deconfined phase is dual to the 5D black-hole thermody- 
namics. The free energy, defined as 

T = E-TS, (2.8) 

is identified with the black-hole on-shell action; as usual, the energy E and entropy S are 
identified with the black-hole mass, and one fourth of the horizon area in Planck units, 
respectively. 

The thermal gas and black-hole solutions with the same temperature differ at 0{r^): 



b{r) = bo{r) 



1+5^ + ... 



,4 



/(0 = l-j^ + --- r^O, (2.9) 



where Q and C are constants with units of energy. As shown in [^] they are related to the 
enthalpy TS and the gluon condensate {TrF'^) : 

TS 22 (TV - (TV 

MSATSyg' ^ 3(47r)2 2A0M^N^ ' ^ ' 

Although they appear as coefficients in the UV expansion, C and G are determined by 
regularity at the black-hole horizon. For T and S the relation is the usual one, 

T = -^, S=^ = AniM^^N^Vs)bHr,). (2.11) 
For G the relation with the horizon quantities is more complicated and cannot be put in a 



simple analytic form. However, as discussed in |1S], for each temperature there exist only 
specific values of Q (each corresponding to a different black hole) such that the horizon is 
regular. 

At any given temperature there can be one or more solutions: the thermal gas is always 
present, and there can be different black holes with the same temperature. The solution 
that dominates the partition function at a certain T is the one with smallest free energy. 



The free energy difference between the black hole and thermal gas was calculated in |19| 
to be: 

^ ^^g_C_ (2.12) 



For a dilaton potential corresponding to a confining theory, like the one we will assume, 



the phase structure is the following |19|: 



*The periods and 3-space volumes of the thermal gas solution are related to the black-hole solution 
values by requiring that the geometry of the two solutions are the same on the (regulated) boundary. See 
lUl for details. 
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1. There exists a minimum temperature r^j„ below which the only solution is the 
thermal gas. 

2. Two branches of black holes ("large" and "small") appear for T > Tmin, but the 
ensemble is still dominated by the confined phase up to a temperature Tc > T^in 

3. At T = Tc there is a first order phase transition to the large black-hole phase. The 
system remains in the black-hole (deconfined) phase for all T > Tc. 

In principle there could be more than two black-hole branches, but this will not happen 
with the specific potential we will use. 



3. Scheme dependence 

There are several sources of scheme dependence in any attempt to solve a QFT. Different 
parametrizations of the coupling constant (here A) give different descriptions. However, 
physical statements must be invariant under such a change. In our case, reparametrizations 
of the coupling constant are equivalent to radial diffeomorphisms as we could use A as the 
radial coordinate. 

In the holographic context, scheme dependence related to coupling redefinitions trans- 
lates into field redefinitions for the bulk fields. As the bulk theory is on-shell, all on-shell 
observables (that are evaluated at the single boundary of space-time) are independent of 
the field redefinitions showing that scheme-independence is expected. Invariance under 
radial reparametrizations of scalar bulk invariants is equivalent to RG invariance. Be- 
cause of renormalization effects, the boundary is typically shifted and in this case field 
redefinitions must be combined with appropriate radial diffeomorphisms that amount to 
RG-transformations. 

Another source of scheme dependence in our setup comes from the choice of the energy 
function. Again we may also consider this as a radial coordinate and therefore it is subject 
to coordinate transformations. A relation between A and E is the /3-function, 

/3(A). (3.1) 



diog^; 



f3 by definition transforms as a vector under A reparametrizations and as a form under E 
reparametrizations. /3(A) can therefore be thought of as a vector field implementing the 
change of coordinates from \ to E and vice-versa. 

Physical quantities should be independent of scheme. They are quantities that are 
fully diffeomorphism invariant. If the gravitational theory had no boundary there would 
be no diffeomorphism invariant quantities, except for possible topological invariants. Since 
we have a boundary, diffeomorphism invariant quantities are defined at the boundary. 

Note that scalar quantities are not invariant. To be invariant they must be scalar and 
constant. We therefore need to construct scalar functions that are invariant under changes 
of radial coordinates. 
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We can fix this reparametrization invariance by picking a very special frame. For 
example choosing the (string) metric in the conformal frame 

ds'^ = e^^ [dr^ + dx^'dx^] , A(r) (3.2) 

or in the domain-wall frame 

ds"^ = du^ + e'^^dx'^dXf, , X{u) (3.3) 

fixes the radial reparametrizations almost completely. In conformal frame, common scalings 
of r,x'^ are allowed, corresponding to constant shifts of A(r). 

Eventually we are led to calculate and compare our results to other ways of calculating 
(like the lattice). Some outputs are easier to compare (for example correlators). Others 
are much harder as they are not invariant (like the value of the coupling at a given energy 
scale) . 

In the UV such questions are well understood. The asymptotic energy scale is fixed 
by comparison to conformal field theory examples. This is possible because the space is 
asymptotically AdSs^. 

The coupling constant is also fixed to leading order from the coupling of the dilaton 
to branes (up to an overall multiplicative factor). Subleading (in perturbation theory) 
redefinitions of the coupling constant and the energy lead to changes in the /3-function 
beyond two loops. 



More in detail, as it has been described in [14, 15 1, the general form of the kinetic term 



for the gauge fields on a brane is expected to be: 

SF2 = e-''Z{R,0Tr[F^] , ? ^ -e'*^ (3.4) 

where Z{R, ^) is an (unknown) function of curvature R and the five-form field strength, ^. 
At weak background fields, Z ~ — ;j + • • • . In the UV regime, expanding near the boundary 



in powers of the coupling A = NcC we obtain, |15] 



Sf2 = N, Tr[F^] J 



(3.5) 



where F{R,S^) is the bulk effective action and R^,^,^ are the boundary values for these 
parameters. Therefore the true 't Hooft coupling of QCD is 

A 



Hooft 



I J Zg(fi*,g*) A ^ 2n 

z(i?„e*)%(i?*,e*)\/e^^ 



(3.6) 



In the IR, more important changes can appear between our A and other definitions as for 
example in lattice calculations. 



"As the dilaton is now not constant there is a non-trivial question: in which frame is the metric AdS. 
In jl^] it was argued that this should be the case in the string frame. The difference of course between the 
string and Einstein frame is subleading in the UV as the coupling constant vanishes logarithmically. But 
this may not be the case in the IR where we have very few criteria to check. In the model we are using we 
impose that the space is asymptotically AdS in the Einstein frame as this is the only choice consistent with 
the whole framework. 
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In the region of strong coupling we know much less in order to be guided concerning 
the correct definition of the energy. We can obtain some hints however by comparing 
with lattice results. ^'^ In particular, based on lattice calculations using the Schrondiger 
functional approach ||5^ , it is argued that at long distance L the 't Hooft coupling constant 
scales as 

Alat ~ e™^ , ~ ^mo++ . (3.7) 

This was based on a specific definition of the coupling constant, and length scale on the 
lattice as well as on numerical data, and some general expectations on the fall-off of corre- 
lations in a massive theory. This suggests an IR f3 function of the form 

^^ = AlogA , A = Aoe™^. (3.8) 

On the other hand our /3-function at strong coupling uses the UV definition of energy, 
log E = Ae (the scale factor in the Einstein frame) , E ^ 1/L and is 



3a- 1 1 
1 + 7 -^ ^ + 

4 a log A 



, A:^ 




(^) 







(3.9) 



where a is a parameter in the IR asymptotics of the potential. The case we consider as 
best fitting YM is a = 2 as then the asymptotic glueball trajectories are linear.. 

Consider now taking as length scale the string scale factor e^" in the IR. Since it 
increases, it is consistent to consider it as a monotonic function of length. From its relation 
to the Einstein scale factor Ag = A^; + | log A and (|3.9D we obtain 

'^^ '^'^ \\og\ + ... . (3.10) 



dAs a — 1 
Therefore if we define as length scale in the IR 

2(7 A 2° 

logL = As L=(e^O"-i (3.11) 

a — 1 

we obtain a running of the coupling compatible with the given lattice scheme. Note however 

2a 

that L = (e^=)''-i cannot be a global choice but should be only valid in the IR. The reason 
is that this function is not globally monotonic. 

We conclude this section by restating that physical observables are independent of 
scheme. But observables like the 't Hooft coupling constant do depend on schemes, and it 
is obvious that our scheme is very different from lattice schemes in the IR. 



^"We would like to thank K. Kajantie for asking the question, suggesting to compare with lattice data, 
and providing the appropriate references. 

^^The string scale factor is not a monotonic function on the whole manifold, and this is the reason 
that it was not taken as a global energy scale. In particular in the UV, e^" decreases until it reaches a 
minimum. The existence of the minimum is crucial for confinement. After this minimum e^" increases and 
diverges at the IR singularity. 
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4. The potential and the parameters of the model 



We will make the following ansatz for the potential^^, 

^(^) = 7? {l + ^oA + FiA^/^ ^i^g (^^ ^ ^ ygA^)] '^'1 , (4.1) 



which interpolates between the two asymptotic behaviors ( |2.2D for small A and (2^) for 
large A, with Q = 2/3 and P = 1/2. Not all the parameters entering this potential have 
physical relevance. Below we will discuss the independent parameters of the model, and 
their physical meaning. 

The normalization of the coupling constant A. As discussed in the previous section, 
the relation between the bulk field A(r) and the physical QCD 't Hooft coupling Xt = 
Qyj^Nc is a priori unknown. In the UV, the identification of the D3-brane coupling to the 
dilaton implies that the relation is linear, and depends on an a priori unknown coefficient 
K, defined as: 

A = K\t. (4.2) 

The coefficient k can in principle be identified by relating the perturbative UV expansion 
of the Yang-Mills /3-function, to the holographic /3-function for the bulk field A: 

22 51 

P{\t) = -/3oA? - /3iA? + . . . /3o = /3i = PI • • • (4.3) 

q q 207 

/3(A) = -6oA2 - 6iA3 + . . . , 6o = ^^o, 6i = 7^1 - .... (4.4) 

o 4 zoo 



The two expressions (O) and ( [4. 41 ) are consistent with a linear relation as in ([4.2|), 



and expanding the identity k/3j(Aj) = /3(KAt) to lowest order leads to: 

K = /3o/5o. (4.5) 

Therefore, to relate the bulk field A to the true coupling At one looks at the linear term in 
the expansion of the potential. More generally, the other /3 function coefficients are related 
by Pn = K"-^^bn, and the combinations 6n/^o^^ = Pu/Pq'^^ are K-independent (however 
they are scheme-dependent for n > 2). 

As discussed in Section |3|, the introduction of the coefficient k amounts to a field 
redefinition and therefore its precise value does not affect physical (scheme-independent) 
quantities. In this sense, k is not a parameter that can be fixed by matching some observable 
computed in the theory. Assuming the validity of the relation (|4.2| ) , we could eventually fix 
K by matching a RG-invariant (but scheme-dependent) quantity, e.g. A at a given energy 
scale. 

However, as we discuss later in this section, rescaling A in the potential (thus changing 
k) affects other parameters in the models, that are defined in the string frame, e.g. the 
fundamental string length ig: if we hold the physical QCD string tension fixed, the ratio 
{is/^) scales with degree —2/3 under a rescaling of k. 



^Further studies of IHQCD with difTerent potentials can be found in |55| 
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An important point to keep in mind, is that the simple hnear relation (O) may be 
modified at strong coupling, but again this does not have any effect on physical observables. 
As long as we compute RG-invariant and scheme-independent quantities, knowledge of the 
exact relationship A = -F(Aj) is unnecessary. 

The AdS scale L This is set by the overall normalization of the potential, and its 
choice is equivalent to fixing the unit of energy. It does not enter dimensionless physical 
quantities. As usual the AdS length at large Nc is much larger than the Planck length 
{£p ~ l/{MpNc^^), independently of the 't Hooft coupling. 

The UV expansion coefficients of ^(A). They can be fixed order by order by matching 
the Yang-Mills /3-function. We impose this matching up to two-loops in the perturbative 
expansion, i.e. O(A^) in /3(A). One could go to higher orders by adding additional powers 
of A inside the logarithm, but since our purpose is not to give an accurate description of 
the theory in the UV, we choose not to introduce extra parameters'^. 

Identifying the energy scale with the Einstein frame scale factor, logE = log6(r), we 
have the relation ( [4. 4] ) between the /3-function coefficients and the expansion parameters 
of y(A), with 

^0 = ^0, vi = Vi^/V2. (4.6) 

The term proportional to V2 in eq. ( [4.1D is needed to reproduce the correct value of the 
quantity 61/69 = l^i/Po = 51/121, which is invariant under rescaling of A. Thus, V2 is not 
a free parameter, but is fixed in terms of Vq and Vi by: 



4 /23 + 366i/6|y 9 61 51 

As explained earlier in this section, when discussing the normalization of the coupling, 
fixing the coefficient Vq is the same as fixing the normalization k through eq. ([4.5|). As we 
argued, the actual value of k should not have any physical consequences, so it is tempting 
to set Vb = 1 by a field redefinition, A — >• A/Vq and eliminate this parameter altogether. 

In fact, most of the quantities we will compute are not sensitive to the value of Vq, but 
for certain quantities, such as the string tension, some extra care is needed. In general, we 
can ask whether two models of the same form ( |2.1| ), but with different potentials V{X) and 
V{X), such that V^(A) = V{aX) for some constant a, lead to different physical predictions. 
As we can change from one model to the other simply by a field redefinition A — )• aX ( this 



has no effect on the other terms in the action in the Einstein frame, eq. (2J) ), clearly 
the two potentials lead to the same result for any physical quantity that can be computed 
unambiguously from the Einstein frame action, e.g. dimensionless ratios between glueball 
masses, critical temperature, latent heat etc. 

However a rescaling of A does affect the string frame metric, since the latter explicitly 
contains factors of A: bs{r) = 6(r)A^/^ 10 thus, under the rescaling A — >• aA, bs{r) — )• 
a^/^bs{r). This means that any dimensionless ratio of two quantities, such that one of 
them remains fixed in the string frame and the other in the Einstein frame, will depend on 



^^Moreover, higher order /3-function coefficients are known to be scheme-dependent. 
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a. An example of this is the ratio is/i, where is is the string length, that we will discuss 
shortly. 

Therefore, we can safely perform a field redefinition and set Vq to a given value, as long 
as we are careful when computing quantities that depend explicitly on the fundamental 
string length. 

Bearing this caveat in mind, we will choose a normalization such that 60 = Po, i-e. 

Vo = 1^0, (4.8) 

so that the normalization of A in the UV matches the physical Yang-Mills coupling. With 
this choice, out of the four free parameters Vi appearing in ( |4.1| ) only Vi and V3 play a 
non-trivial role (V2 being fixed by eq. ( [4.7D ). 

The 5D Planck scale Mp. Mp appears in the overall normalization of the 5D action 
( |2.1| ). Therefore it enters the overall scale of quantities derived by evaluating the on-shell 
action, e.g. the free energy and the black-hole mass. It also sets the conversion factor 
between the entropy and horizon area. Mp cannot be fixed directly as we lack a detailed 
underlining string theory for YM. To obtain quantitative predictions, Mp must be fixed in 
terms of the other dimension-full quantity of the model, namely the AdS scale 1. As shown 



in |19| this can be done by imposing that the high-temperature limit of the black-hole free 
energy be that of a free gluon gas with the correct number of degrees of freedom^^. This 
requires: 

{Mpif = (4.9) 

The string length. In the non-critical approach the relation between the string length 
is and the 5D Planck length (or the AdS length i) is not known from first principles. 



The string length does not appear explicitly in the 2-derivative action (2J), but it enters 
quantities like the static quark-antiquark potential. The ratio is/i can be fixed phenomeno- 
logically to match the lattice results for the confining string tension. 

More in detail, the relation between the fundamental and the confining string tensions 
Tf and a is given by: 

a = Tfb^{n)X^/^{n), (4.10) 

where r^, is the point where the string frame scale factor, bs{r) = 6(r)A^/^(r), has its 
minimum. Fixing the confining string tension by comparison with the lattice result we 
can find Tf (more precisely, the dimensionless quantity Tfi'^, since the overall scale of the 
metric depends on i). The string length is in turn given by is/i = ^lixTji? . 



As is clear from eq. ( 4.10 ), rescaling A — t- a A, keeping the value of the QCD string 
tension a and of the AdS scale I fixed, affects the fundamental string length in AdS units 

as Isji — )• a~'^^^{is/ i)- Therefore two models a and h, defined in the Einstein frame by 



^*Note that this is conceptually different from the A/" = 4 case. There, near the boundary, the theory 
is strongly coupled and this number must be calculated in string theory. It is different by a factor of 3/4 
from the free sYM answer. Ifere near the boundary the theory is free. Therefore the number of degrees of 
freedom can be directly inferred. 



- 15 - 



eq. (|2.1[) , but with potentials related by Vf,(A) = Va{aX), must have different fundamental 
string tensions in order to reproduce the same result for the QCD string tension. The 
quantity is/(^ therefore depends on the value of Vq. 



Integration constants. Besides the parameters appearing directly in the gravitational 
action, there are also other physically relevant quantities that label different solutions to 
the 5-th order system of field equations. Any solution is characterized by a scale A, the 
temperature T and a value for the gluon condensate Q, that correspond to three of the five 
independent integration constants. 

Regularity at the horizon fixes ^ as a function of T, so that effectively the gluon 
condensate is a temperature-dependent quantity. 

The quantity A controls the asymptotic form of the solution, as it enters the dilaton 
running in the UV: A ~ — (60 log rA)~^. It can be defined in a reparametrization invariant 
way as: 





1 


exp 





A = rMimU(A) tu^y (4-11) 

and it is fixed once we specify the value of the scale factor 6(A) at a given Aq- 

Every choice of A corresponds to an inequivalent class of solutions, that differ by UV 
boundary conditions. Each class is thermodynamically isolated, since solutions with dif- 
ferent A's have infinite action difference. Thus, in the canonical partition sum we need to 
consider only solutions with a fixed value of A. However, this choice is merely a choice of 
scale, as solutions with different A's will give the same predictions for any dimensionless 
quantity. In short, A is the holographic dual to the QCD strong coupling scale: it is defined 
by the initial condition to the holographic RG equations, and does not affect dimensionless 
quantities such as mass ratios, etc. Therefore, as long as all solutions we consider obey the 
same UV asymptotics, the actual value of A is immaterial, since the physical units of the 
system can always be set by fixing i. 

To summarize, the only nontrivial phenomenological parameters we have at our disposal 
are Vi and V3 appearing in ( |4.1| ). The other quantities that enter our model are either 
fixed by the arguments presented in this section, or they only affect trivially (e.g. by 
overall rescaling that can be absorbed in the definition of the fundamental string scale) the 
physical quantities. 

In the next section we present a numerical analysis of the solutions and thermody- 
namics of the model defined by eq. (|4.l|), and show that for an appropriate choice of the 
parameters it reproduces the lattice results for the Yang-Mills deconfinement transition 
and high-temperature phase as well as the zero temperature glueball data. 
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The remaining two are the value /(O) which should be set to one for the solution (2.7) to obey the right 



UV asymptotics, and an unphysical degree of freedom in the reparametrization of the radial coordinate. 
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5. Matching the thermodynamics of large- Ai'c YM 



Assuming a potential of the form ( [4.1D , we look for values of the parameters such that the 
thermodynamics of the 5D model match the lattice results for the thermodynamics of 4D 
YM. As explained in Section |3|, we set Vq and V2 as in eqs. ( |4.8| ) and (^]^), respectively, 
with 60 = /So = 22/3(47r)-2. 

We then vary Vi and V3 only. We fix these parameters by looking at thermodynamic 
quantities corresponding to the latent heat per unit volume, and the pressure at one value 
of the temperature above the transition, which we take as 2Tc. 

It is worth remarking that Vi and V3 are phenomenological parameters that we use 
to fit dimensionless QCD quantities. The single (dimension-full) parameter of pure Yang- 
Mills, the strong coupling scale, is an extra input that fixes the overall energy scale of our 
solution. 

Using the numerical method explained in [22|, for each set of parameters (Vi, V3) we 
numerically generate black-hole solutions for a range of values of \h, then from the metric 
at the horizon and its derivative we extract the temperature and entropy functions T{Xh) 
and S{Xh), and the function J-{Xh) from the integrated form of the first law. 



dXhS{Xh 



dTiXh 



(5.1) 



Here S{Xh) is given by ( |2.11| ) and both the large black hole and small black-hole branches 
are needed in order the get the full result for the free energy. This is because the integral 
in (5.1) extends to +00, entering deeply in the small black-hole branch. 

The behavior of the thermodynamic functions is shown in Figures ||, ^ and ^, for the 
best fit parameter values that we discuss below. One can see the existence of a minimal 
temperature Tmin = T{Xmin), and a critical value Ac where T changes sign. The resulting 
function J-{T) is shown in Figure |l[ 




Figure 1: The Free energy density (in units of Tc) as a function of T/Tc, for Vi = 14 and 
V3 = 170. The vertical lines correspond to the critical temperature (solid) and the minimum 
black-hole temperature (dashed). 

The phase transition is first order, and the latent heat per unit volume Lh, normalized 
by N^T^, is given by the derivative of the curve in Fig. || at T/Tc = 1. Equivalently, Lh 
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is proportional to the jump in the entropy density s = S/V3 at the phase transition from 
the thermal gas (whose entropy is of 0(1), in the limit Nc ^ 00) to the black hole (whose 
entropy scales as in the same limit): thus, in the large Nc limit, 

Lh = TAs ~ T,s{Xc) (5.2) 

up to terms of 0(1 /N^). 

To fix Vi and V3 we compare our results to the data of G. Boyd et al. |5^. The relevant 
quantities to compare are the dimensionless ratios p(T)/T^, e(T)/T^ and s(T)/T^, where 
p = J- jVj, is the pressure, and e = p + Ts is the energy density. Lattice results for these 
functions are available in the range T = Tc ^ 5Tc, and can be seen in Figure 7 of |5^]. The 
analysis of ||5^ correspond to Nc = 3, but one expects that the thermodynamic functions 
do not change to much for large Nc^^. 

T 



1.4 - 




Figure 2: Temperature in units of Tmin, as a function of Xh, for Vi = 14 and V3 — 170. The 
dashed horizontal and vertical lines indicate the critical temperature and the critical value of the 
dilaton field at at the horizon 

An additional quantity of relevance is the value for the "dimensionless" latent heat per 
unit volume, L/^/T^ which for large Nc was found in H to be {Lh/T^)iat = O.SliV^^ The 
resuh for N^ = 3 is slightly lower ( ~ 0.28N^). 

As already noted in |l9[ , the qualitative features of the thermodynamic functions 
are generically reproduced in our setup: the curves 3p{T)/T^, e{T)/T'^ and 3s(T)/4r^ 
increase starting at Tc, then (very slowly) approach the constant free field value tt'^N^/15 
(given by the Stefan-Boltzmann law) as T increases. By computing the thermodynamic 
functions for various sets of values of Vi and V3 we obtain that: 

1) Vi roughly controls the height reached by the curves p{T)/T'^, e{T)/T'^ and s{T)/T'^ 
at large T/Tc (~ a few): for larger Vi the curves approach the free field limit faster; 

2) V3 does not affect much the height of the curves at large T/Tc, but on the other 
hand it changes the latent heat, which is increasing as V3 decreases. 

The best fit corresponds to the values 

^®See e.g. |^^, in which results for A'^c = 8 do not different significantly from those for A^c = 3 as well as 
the recent high-precision data by Panero, ^A. 
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F 




Figure 4: Entropy density in units of Tmin, as a function of Xh 



Vi = 14 Vi = 170. (5.3) 

Below we discuss the values of various physical quantities ( both related to thermodynamics, 
and to zero-temperature properties) obtained with this choice of parameters. 

5.1 Latent heat and equation of state 

The comparison between the curves p{T) /T'^, e{T)/T'^ and s(T) /T^ obtained in our models 
with ( |5.3D , and the lattice results ||5^ is shown in Figure ||. The match is remarkably good 
for Tc < T < 2Tc, and deviates slightly from the lattice data in the range up to STc. 
The latent heat we obtain is: 

Lh/T^ = 0.31Nl (5.4) 

which matches the lattice result for Nc ^ oo |5^. 

An interesting quantity is the trace anomaly (e — 3p)/T^, (also known as interaction 
measure), that indicates the deviation from conformality, and it is proportional to the 
thermal gluon condensate. The trace anomaly in our setup is shown, together with the 
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Figure 5: Temperature dependence of the dimensionless thermodynamic densities s/T^ (hght 
blue), p/T^ (dark blue) and e/T** (green), normalized such that they reach the common limiting 
value 7r^/15 (dashed horizontal line) as T — )■ oo. The dots correspond to the lattice data for iVc = 3 
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Figure 6: The trace anomaly as a function of temperature in the deconfined phase of the holo- 
graphic model (solid line) and the corresponding lattice data [|6j for Nc — i (dots). The peak in 
the lattice data slightly above Tc is expected to be an artifact of the finite lattice volume. In the 
infinite volume limit the maximum value of the curve is at Tc, and it equals Lh/N^T^- 



corresponding lattice data, in figure (|6|), and the agreement is again very good. Our results 
agree even better with recent high-precision lattice calculations of the thermodynamics 
functions done by Panero at different values of Nc up to Nc = 8, |5^. In figure ^ a 
comparison (taken from [^]) of the normalized interaction measure with lattice results for 
different Nc is shown. 

We also compute the specific heat per unit volume Cy, and the speed of sound Cg in 
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Figure 7: The rescaled trace anomaly (so that it is iVc-independent) as a function of temperature 
in the deconfined phase of the holographic model (solid line) and the corresponding recent high 



precision lattice data taken from 58 for different N^- The errors shown are statistical only. 

C, 




Figure 8: The specific heat (divided by T^^), as a function of temperature, in the deconfined phase 
of the holographic model. 



the deconfined phase, by the relations 



(5.5) 



These are shown in Figures |8| and ^ respectively. The speed of sound is shown together 
with the lattice data, and the agreement is remarkable. 

5.2 Glueball spectrum 



In 1 14], the single phenomenological parameters of the potential was fixed by looking at the 
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Figure 9: The speed of sound in the deconfined phase, as a function of temperature, for the 
holographic model (solid line) and the corresponding lattice data for Nc = 3 (dots). The 
dashed horizontal line indicates the conformal limit — 1/3. 



zero-temperature spectrum, i.e. by computing various glueball mass ratios and comparing 
them to the corresponding lattice results. The masses are computed by deriving the ef- 
fective action for the quadratic fluctuations around the background, and subsequently 
reducing the dynamics to four dimensions. 

The associated thermodynamics for this potential was studied in which was in 
qualitative agreement with lattice QCD results, but not in full quantitative agreement. This 
is due to the fact that the thermodynamics depends more on the details of the potential 
than the glueball spectrum for the main Regge trajectories. Here we use the potential 



(4.1), but with the two phenomenological parameters Vi and V3 already determined by the 
thermodynamics (5.3). 

The glueball spectrum is obtained holographically as the spectrum of normalizable 
fluctuations around the zero-temperature background. As explained in the introduction, 
and motivated in here we consider explicitly the 5D metric, one scalar fleld (the 

dilaton), and one pseudoscalar fleld (the axion). As a consequence, the only normalizable 
fluctuations above the vacuum correspond to spin and spin 2 glueballs^^ (more precisely, 
states with J^'^ = 0++, 0"+, 2++), each species containing an inflnite discrete tower of 
excited states. 

In 4D YM there are many more operators generating glueballs, corresponding to dif- 
ferent values of J^^ , that are not considered here. These are expected to correspond 
holographically to other fields in the noncritical string spectrum (e.g. form fields, which 
may yield spin 1 and CP-odd spin 2 states) and to higher string states that provide higher- 
spin glueballs. As the main focus is in reproducing the YM thermodynamics in detail 
rather than the entire glueball spectrum, we choose not to include these states^^. There- 
fore we only compare the mass spectrum obtained in our model to the lattice results for 

^^Spin 1 excitations of the metric can be shown to be non-normalizable. 

further reason is that, unlike the scalar and (to some extent) the pseudoscalar sector that we are 
considering, the action governing the higher Regge slopes is less and less universal as one goes to higher 
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the lowest 0^^,0 ^,2^+ glueballs and their available excited states. These are limited to 
one for each spin species, and none for the spin 2, in the study of |61|, which is the one 
we use for our comparison. This provides two mass ratios in the CP-even sector and two 
in the CP-odd sector. 

The glueball masses are computed by first solving numerically the zero-temperature 
Einstein's equations, by setting /(r) = 1, and using the resulting metric and dilaton to 
setup an analogous Schrodinger problem for the fluctuations, [14|. The results for the 
parity-conserving sector are shown in Table ||, and are in good agreement with those 
reported by |61] for Nc = 3, whereas the results reported by |62] for large Nc are somewhat 
larger. The CP-violating sector (axial glueballs) will be discussed separately. 



Table 1: Glueball Masses 





HQCD 


A^c = 3 


iVc = oo H] 


mQ*++ /mQ++ 


1.61 


1.56(11) 


1.90(17) 


1712++ /mQ++ 


1.36 


1.40(4) 


1.46(11) 



We should add that there are other lattice studies (see e.g. |Q) that report additional 
excited states. Our mass ratios offer a somewhat worse fit of the mass ratios found in |63| 
(whose results are not entirely compatible with those of for the states the two studies 
have in common). We should stress however that reproducing the detailed glueball spec- 
trum is secondary here since the main focus is thermodynamics. However, the comparison 
of our spectrum to the existing lattice results shows that our model provides a good global 
fit to 4D YM also with respect to quantities beyond thermodynamics. 

Unlike the various mass ratios, the value of any given mass in AdS-length units (e.g. 
771-0++^) does depend on the choice of integration constants in the UV, i.e. on the value of 
bjjv and Xuv- Therefore its numerical value does not have an intrinsic meaning. However 
it can be used as a benchmark against which all other dimension-full quantities can be 
measured (provided one always uses the same UV b.c. ). On the other hand, given a fixed 
set of initial conditions, asking that mo++ matches the physical value (in MeV) obtained 
on the lattice, fixes the value of i hence the energy unit. 



5.3 Critical temperature 

The thermodynamic quantities we have discussed so far, are dimensionless ratios, in units 
of the critical temperature. To compute Tc, we need an extra dimension- full quantity which 
can be used independently to set the unit of energy. In lattice studies this is typically the 
confining string tension a in the T = vacuum, with a value of around (440Mey)^, and 
results are given in terms of the dimensionless ratio Tc/^/a. In our case we cannot compute 
a directly, since it depends on the fundamental string tension, which is a priori unknown. 
Instead, we take the mass ?Tio++ of the lowest-lying glueball state as a reference. 

masses. Only a precise knowledge of the underline string theory is expected to provide detailed information 
for such states. 
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We compute mo++ with the potential (|4.1|) , with Vi and V3 fixed as in ( p^) , then 
compare Tc/mo++ to the same quantity obtained on the lattice. For the lattice result, we 
take the large Nc result of Tc/y/a = 0.5970(38), and combine it with the large Nc 
result for the lowest-lying glueball mass [^, mQ^^/^/a = 3.37(15). The two results are in 
fair agreement, without need to adjust any extra parameter: 



0.167, — = 0.177(7) . (5.6) 



—] =0.167, , 

^oJhQCD \'^0 J lattice 

In physical units, the critical temperature we obtain is given by 

Tc = 0.56 ^/a = 247 MeV. (5.7) 

5.4 String tension 

The fundamental string tension Tf = cannot be computed from first principles in our 
model, but can be obtained using as extra input the lattice value of the confining string 
tension cr, at T = 0. The fundamental and confining string tensions are related by eq. 



As for the critical temperature, we can relate Tf to the value of the lowest-lying glueball 
mass, by using the lattice relation y/a = -^37^ |62]. Since what we actually compute 



numerically is mo++^, this allows us to obtain the string tension Tf (and fundamental 
string length ig = Ij ^ ^vrTj in AdS units: 

Tfl^ = 0.19, is/i = 0.15 . (5.8) 

This shows that the fundamental string length in our model is about an order of magni- 
tude smaller than the AdS length. The meaning of this fact is a little more complicated 
conceptually, as the discussion in |15] indicates. Also, we should stress that, as discussed 



in Section ^ this result depends on our choice of the overall normalization of A: changing 
the potential by A — )• kA will yield different numerical values in ( |5.8D without affecting the 
other physical quantities. 

Another related observable is the spatial string tension. It is calculated from the 
expectation value of the rectangular Wilson loop which stretches in spatial dimensions 
only. This has been calculated on the lattice [^], as well as using the high-temperature 
(resumed) perturbative expansion plus a zero-temperature calculation of the string tension 



in three-dimensional YM theory, |65]. The two calculation agree reasonably well. 

The spatial string tension at finite temperature can be calculated in IHQCD, [^] by 
calculating the relevant Wilson loop. Very good agreement was found with the lattice 
calculations, especially at temperatures not far from the phase transition. 

Finally, several calculations of quark- antiquark potentials exist. At zero temperature 
the long distance asymptotics of the quark potential was calculated in [0] and used to 
classify the dilaton potentials as a function of the confinement property. The full quark 



potential including the short distance behavior was computed in [S7|. There a comparison 



to the Cornell potential was done as well as with quarkonium spectra finding excellent 
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agreement with data. The issue of quarkonium potentials from IHQCD-hke theories was 
also recently discussed in [|6^ ]. 

Finally the Polyakov loop was recently computed ||6^ in similar Einstein dilaton models 
that were studied first in [p^]. 

5.5 CP-odd sector 

The CP-odd sector of pure Yang-Mills is described holographically by the addition of a 
bulk pseudoscalar field a(r) (the axion) with action:^^ 



±v±p 



d^xZ{X)^{df'a){d^a) . (5.9) 



The field a(r) is dual to the topological density operator TrFF. The prefactor Z{\) is a 
dilaton-dependent normalization. The axion action is suppressed by a factor with 
respect to the action ( |2.lD for the dilaton and the metric, meaning that in the large-A^c 
limit one can neglect the back-reaction of the axion on the background. 



As shown in |14], requiring the correct scaling of a(r) in the UV, and phenomenologi- 



cally consistent axial glueball masses, constrain the asymptotics of Z(X) as follows: 

Z{X) ~ Zo , A ^ 0; Z(A) ~ A ^ oo, (5.10) 

where Zq is a constant. As a simple interpolating function between these large- and small- 
A asymptotics we can take the following: 

ZiX) = Zo{l+CaX^). (5.11) 

The parameter Zq can be fixed by matching the topological susceptibility of pure Yang- 
Mills theory, whereas can be fixed by looking at the axial glueball mass spectrum. 

Axial glueballs. As in 1 14] , we can fix Ca by matching to the lattice results the mass ratio 

mo ^/mo++ between the lowest-lying axial and scalar glueball states. This is independent 

of the overall coefficient Zq in ( 5.11] ). The lattice value uiq ^/mo++ = 1.49 |61] is obtained 



for: 

Ca = 0.26. (5.12) 

With this choice, the mass of the first excited axial glueball state is in good agreement 
with the corresponding lattice result []6l| : 

^ -2.10 (a^:^) =2.12(10). (5.13) 



"^0++ / hQCD V "10++ / lattice 



^^This action was justified in p^ , [l^ . The dilaton dependent coefficient Z{X) is encoding both the 
dilaton dependence as well as the UV curvature dependence of the axion kinetic terms in the associated 
string theory. We cannot determine it directly from the string theory, but we pin it down by a combination 
of first principles and lattice input, as we explain further below. 
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Topological Susceptibility. In pure Yang-Mills, the topological x susceptibility is de- 
fined by: 



(5.14) 



where E{0) is the vacuum energy density in presence of a ^-parameter. E[6) can be 
computed holographically by solving for the axion profile a(r) on a given background, and 
evaluating the action (^) on-shell. 

In the deconfined phase, the axion profile is trivial, implying a vanishing topological 
susceptibility [^]. This is in agreement with large- iVc arguments and lattice results pT| . 

In the low-temperature phase, the axion acquires a non-trivial profile. 



a{r) = auv 



F{r) 



Fir) = 



dr 



Z{X{r))e^Mr) ' 



(5.15) 



This profile is shown, for the case at hand, in Figure |l^, where the axion is normalized to 
its UV value. 




100 200 300 400 500 600 



Figure 10: Axion profile in the radial direction. The a;-axis is taken to be the energy scale, 
E{r) — Eob{r), where the unit Eq is fixed to match the lowest glueball mass. 



The topological susceptibility is given by p4[| : 



dr 



-M(r)z{r) 



(5.16) 



where Z{r) = Z(A(r)). Evaluating this expression numerically with Z{X) as in ( 5.11| ), and 
Ca = 0.26 (to match the axial glueball spectrum ), we can determine the coefficient Zq 
by looking at the lattice result for x- Nc = 3, [^0| obtained x = (WlMeV)^ , which 
requires Zq = 133. 

In Table Q we present a summary of the various physical quantities discussed in this 
section, as obtained in our holographic model, and their comparison with the lattice results 
for large Nc (when available) and for Nc = 3. The quantities shown in the upper half of the 
table are the ones that were used to fix the free parameters (reported in the last column) 
of the holographic model. 
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5.6 Coupling normalization 



Finally, we can relate the field A(r) to the running 't Hooft coupling. All other quantities 
we have discussed so far are scheme-independent and RG-invariant. This is not the case 
for the identification of the physical YM 't Hooft coupling, which is scheme dependent. 
In the black- hole phase we can take = X{rh) as a measure of the temperature- 



dependent coupling. In figure 11 we show A/^ as a function of the temperature in the range 
Tc to 5Tc. 




Figure 11: The coupling at the horizon as a function of temperature in the range T^-BTc- 

As a reference, we may take the result of @, that found g^{5Tc) ~ 1.5 for Nc = 3, 
which translates to Xt{5Tc) — 5. On the other hand, if we make the assumption that 
the identification A = Aj is valid at all scales (not only in the UV), we find in our model 
At(5Tc) ~ 0.04 (see Figure 0), i.e. a factor of 100 smaller than the lattice result. 

This discrepancy is almost certainly due to the identification ( [4. 2]) being very different 
from lattice at strong coupling. 



6. Bulk viscosity 



The bulk viscosity C is an important probe of the quark-gluon plasma. Its profile as 
a function of T reveals information regarding the dynamics of the phase transition. In 
particular, both from the low-energy theorems and lattice studies |]29| , 30, there is 
evidence that C increases near Tc. 

For a viscous fluid the shear r] and bulk C viscosities are defined via the rate of entropy 
production as 



ds 
di 



1 
T 



diVj + djVi 



(6.1) 



Therefore, in a holographic setup, the bulk viscosity can be defined as the response of 
the diagonal spatial components of the stress-energy tensor to a small fluctuation of the 
metric. It can be directly related to the retarded Green's function of the stress-energy 
tensor by Kubo's linear response theory: 



C = lim -ImGRiu, 0), 



(6.2) 
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Table 2: Collected in this table is the complete set of physical quantities that we computed in 
our model and compared with data. The upper half of the table contains the quantities that we 
used as input (shown in boldface) for the holographic QCD model (HQCD). Each quantity can be 
roughly associated to one parameter of the model (last column) . The lower half of the table contains 
our "postdictions" (i.e. quantities that we computed after all the parameters were fixed) and the 
comparison with the corresponding lattice results. The value we find for the critical temperature 
corresponds to = 24:7 MeV. 







lattice iVc — o 


lattice N[. — )■ oo 


Parameter 


mQ++ 1 y/d 
mQ-+ /mQ++ 

X 


1.2 
0.31 

7rV45 
6.6 1 
1.49 

(191Mey)^ 


1.2 

0.28 156] 

^2/45 
o.oo |bl| 
1.49 H 


0.31 [59[ 

7rV45 
6.6i p2| 


VI = 14 
V2, = 170 

Mpl = [457r2]-i/3 

ts/t = 0.15 
Ca = 0.26 
Zo = 133 


Tc/mQ++ 


0.167 




0.177(7) 




m2++ / mQ++ 


1.61 
1.36 


1.56(11) 
1.40(4) 


1.90(17) 
1.46(11) 




mo*-+/mo++ 


2.10 


2.12(10) 







where Gr{w,p) is the Fourier transform of retarded Green's function of the stress-energy 
tensor: 

3 



Gr{w,p) 



(6.3) 



A direct computation of the RHS on the lattice is non-trivial as it requires analytic contin- 
uation to Lorentzian space-time. In refs. |29],|3C] the low energy theorems of QCD, as well 
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as (equilibrium) lattice data at finite temperature were used in order to evaluate a partic- 
ular moment of the spectral density of the relevant correlator, using a parametrization of 
the spectral density via two time-dependent constants, one of which is the bulk viscosity a 
relation for their product was obtained as a function of temperature. This can be converted 
to a relation for assuming the other constant varies slowly with temperature. 

The conclusion was that (^/s increases near Tc- Another conclusion is that the fermionic 
contributions to C ^-re small compared to the glue contributions. 



The weak point of the approach of |3C], is that it requires an ansatz on the spectrum 
of energy fluctuations, and further assumptions on the other parameters, which are not 
derived from first principles. 



A direct lattice study of the bulk viscosity was also made in fS^]. Here, the result is 



also qualitatively similar 12. However, the systematic errors in this computation are large 
especially near Tc, mostly due to the analytic continuation that one has to perform after 
computing the Euclidean correlator on the lattice. 



The results of references , |Q and the assumptions of the lattice calculation have 



been recently challenged in |71] 



6.1 The holographic computation 

The holographic approach offers a new way of computing the bulk viscosity. In the holo- 
graphic set-up, C is obtained from ( |6.2D . Using the standard AdS/CFT prescription, the 
two point-function of the energy-momentum tensor can be read off from the asymptotic 
behavior of the metric perturbations Sg^i, . This is similar in spirit to the holographic com- 
putation of the shear viscosity fT^ , but it is technically more involved. A recent treatment 
of the fluctuation equation governing the scalar mode of a general Einstein-Dilaton system 



can be found in [^]. Here, we follow the method proposed by [74|. 

As explained in one only needs to examine the equations of motion in the gauge 
r = <I>, where the radial coordinate is equal to the dilaton. In our type of metrics, the 
applicability of this method requires some clarifications, that we provide in [^^. Using 
SO (3) invariance and the five remaining gauge degrees of freedom the metric perturbations 
can be diagonalized as 

Sg = diag{goo, gii, gu, gu, 955), (6.4) 

where 

goo = f[l + hoo{ip)e-''% gn = e^^[l + hu{'p)e-''% , (6.5) 

T 

where the functions A and B emerge from the metric 

ds^ = e^M^X-fdt" + dxmdx^) + e2^(^)^. (6.6) 

Here, the fluctuations are taken to be harmonic functions of t while having an arbitrary 
dependence on (p. 



555 = V[l + ^55(v=)e-*'^*] 
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The bulk viscosity depends only on the correlator of the diagonal components of the 
metric and so it suffices to look for the asymptotics of hn. Interestingly, in the r = $ 
gauge this decouples from the other components of the metric and satisfies the following 
equation^'^ 

- 7) - (-^"' + 97^ - ^) = » • f-^' 

One needs to impose two boundary conditions. First, we require that only the infalling 
condition survives at the horizon: 

hu ^ Cb{(ph- , ^p^(ph, (6.8) 

where Cb is a normalization factor. The second boundary condition is that /in has unit 
normalization on the boundary: 

/ill 1, 93 -00. (6.9) 



Having solved for /in (99), Kubo's formula (6.2) and a wise use of the AdS/CFT prescription 
to compute the stress-energy correlation function [ 74 1 determines the ratio of bulk viscosity 
as follows. 

The AdS/CFT prescription relates the imaginary part of the retarded Ta Green's 
function to the number flux of the /in gravitons J- [74]: 

ImGRiu;,0) = --f— (6.10) 

where the flux can be calculated as the Noether current associated to the U{1) symmetry 
/ill — ^ e*^/iii in the gravitational action for fluctuations. One finds, 

^ = i^jfT^lKik'n - hnh*n]. (6.11) 

As is independent of the radial variable, one can compute it at any (/?, most easily near 
the horizon, where /in takes the form ( |6.8[ ). Using also the fact that {dA/dip){iph) = 
-8Vi^h)/9V'i^h), one finds 

-^(0) = |0|c.(0)Pe^-^(-)^^. (6.12) 

Then, ( |6.2D and (6.10) determine the ratio of bulk viscosity and the entropy density as. 

In the derivation we use the Bekenstein-Hawking formula for the entropy density, s = 
exp3A(v?fe)/4G5. 

To find C we need to find only in the limit cj — t- 0. The computation is performed by 
numerically solving equation (^^) with the appropriate boundary conditions. There are 
two separate methods that one can employ to determine the quantity 0^: 

^"DifTerence in the various numerical factors in this equation w.r.t ||7^ is due to our difTerent normahzation 
of the dilaton kinetic term. 
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Figure 12: Plot of C/s (continuous line) calculated in Improved Holographic QCD model. This 
is compared with the lattice data of that are shown as boxes. The horizontal dashed line is 
indicating the (universal) value of j for comparison. 



1. One can solve ( |6.7| ) numerically with a fixed ijj/T, but small enough so that Cf, reaches 
a fixed value. The method is valid also for finite values of cj. From a practical point 
of view, it is easier to solve ( |6.7| ) with the boundary condition ( |6.8| ) with a unit 
normalization factor, read off the value on the boundary /in (—00) from the solution 
and finally use the symmetry of ( |6.7| ) under constant scalings of hu to determine 
\ch\ = 1/1 /ill (-00) I . 



2. An alternative method of computation that directly extracts the information at = 



follows from the following trick Instead of solving (6^) for small but finite w, 
one can instead solve it for uj = This is a simpler equation, yet complicated 
enough to still evade analytic solution. Let us call this solution h\i. One numerically 
solves it by fixing the boundary conditions on the boundary: /i]'i(— 00) = 1 and the 
derivative dh^i/d(p{— 00) is chosen such that hu is regular at the horizon. Matching 
this solution to the expansion of ( |6.8| ) for small oj than yields |cb| = h\i{iph). 

Both methods were used to obtain C,/ s as a function of T and checked that they yield 
the same result. As explained in most of the thermodynamic observables are easily 
computed using the method of scalar variables, [76|. 



The results are presented in figure 12. This figure gives a comparison of the curve 
obtained by the holographic calculation sketched above by solving ( |6.7D and the lattice 
data of ||3^. We also show r//s = l/47r in this figure for comparison. The result is 
qualitatively similar to the lattice result where C,/ s increases as T approaches Tc, however 
the rate of increase is slower than the lattice. As a result, we obtain a value C/s{Tc) ~ 0.06 



that is an order of magnitude smaller than the lattice result|32| which is 0.8. Note however 
that the error bars in the lattice evaluation are large near Tc and do not include all possible 
systematic errors from the analytic continuation. 
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We should note the fact that the holographic calculation gives a smaller value for the 
bulk viscosity near Tc than the lattice calculation is generic and has been found for other 
potentials with similar IR asymptotics, [74|. The fact that the value of C/s near Tc is 
correlated with the IR asymptotics of the potential will be shown further below. 

Another fact that one observes from figure [T^ is that (/s vanishes in the high T limit. 
This reflects the conformal invariance in the UV and can be shown analytically as follows. 
C/s is determined by formula ( |6.13| ). In the high T limit, (corresponding to Xh — )• 0, 
near the boundary), the fluctuation coefficient |cf,| — )■ 1. This is because of the boundary 
condition /in (A = 0) = 1. We use the relation between T and Xh in the high T limit fH 



A;,^ (/3olog(7rr/A))-i. (6.14) 
Substitution in ( |6.13| ) leads to the result, 

c 



s 



2 1 , , 

— )• — ^ , as T ^ DO. (6.15) 

large 277r log^ (vrT/A) ' ^ ' 



As s itself diverges as in this limit - it corresponds to an ideal gas - we learn that C, 
also diverges as r^/log^(T). Divergence at high T is expected from the bulk- viscosity of 
an ideal gas. We do not expect however the details of the asymptotic result to match with 
the pQCD result, for the same reasons that the shear-viscosity-to-entropy ratio does not, 
flSl] . We note however, that the asymptotic T-dependence is very similar to the pQCD 



result, [[75| : 



C/s oc log-^(7rr/A) log-^ log(7rr/A) . (6.16) 

6.2 Holographic explanation for the rise of (/s near and the small black- hole 
branch 

With the same numerical methods, one can also compute the ratio C/s on the small black- 
hole branch. As this solution has a smaller value of the action than the large black-hole 
solution, it is a subleading saddle point in the phase space of the theory, hence bears no 
direct significance for an holographic investigation of the quark- gluon plasma. However, 
as we show below, the existence of this branch provides a holographic explanation for the 
peak in ^/s in the quark-gluon plasma, near Tc. 

From the practical point of view, we find the second numerical method above (solving 
the fiuctuation equation at u = 0) easier in the range of Xh that corresponds to the small 
black hole. The result is shown in figure ^ (a). The presence of two branches for T > Tmin, 
is made clear in this figure. See also fig[l^ (b) for the respective ranges of Xh that correspond 



to small and large BHs. In fig |13| (a), ("/s on the large BH branch is depicted with a solid 
curve and the small BH branch is depicted with a dashed curve. We observe that (^/s keeps 
increasing on the large-BH branch as T is lowered, up to the temperature Tmin where the 
small and large BH branches merge^^. On the other hand, on the small BH branch (/s 
keeps increasing as the T is increased, up to a certain T^ax that lies between T^m and Tc, 



see figure 14. From this point onwards, C/s decreases with increasing T. 



^^As far as the thermodynamics of the gluon plasma is concerned, the temperatures below Tc (on the 
large BH branch) has little importance, because for T <Tc the plasma is in the confined phase. 
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(a) (b) 

Figure 13: (a) Numerical evaluation of C/?? both on the large-BH branch (the solid curve) and on 
the small BH branch (the dashed curve). T„i denotes Tmin- (b) The two branches of black- hole 
solutions, that correspond to different ranges of Xh- The large BH corresponds to Xh < Xmin and 
the small BH corresponds to X^ > Xmm- 



Tc 



Tnjax 



Trhin 



0.95 0.96 0.97 0.98 0.99 1.00 1.01 1.02 



T/Tc 



Figure 14: An inset from the figure [TbI around the maximum of (/s. 

A simple fact that can be proved analytically is that the derivative of C/s diverges at 
Tmin- This is also clear from figure p!^. This is shown by inspecting equation ( |6.13| ). The 
T derivative is determined as d/dT = {dT / dXh)d/ dX^. Whereas the derivative w.r.t X^ is 
everywhere smooth^^, the factor dT/dXh diverges at Tmin by definition, see figure [l^ (b). 

Therefore, the presence of a Tmin where the large and the small black holes meet, in 
other words, the presence of a small black-hole branch is responsible for the increase of C/s 
near Tmin- As in most of the holographic constructions that we analyzed, and specifically 
in the example we present here, T^ and Tmin are close to one another, this fact implies a 
rise in the bulk viscosity near T^. This proposal, combined with the fact that the existence 
of a small BH branch and color confinement in the dual gauge theory at zero T are in one- 



22 



Note that Cb is also a function of A^. As both the fluctuation equation (6.7) and the boundary conditions 



are smooth at Xh ~ Xmin, one concludes that Cb also is smooth at this point. 
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to-one correspondence []19| , suggests that in confining large-N gauge theories, there will be 
a peak in the ratio (/s close to Tc- 

Another fact that can be shown analyticahy is that C,/ s asymptotes to a finite value 
as T — 7- oo in the smah black- hole branch We find that, 



1 

^ — , as oo. (6.17) 

small 



As the entropy density vanishes in this limit [19|, we conclude that C, should vanish with 
the same rate. 

For a general potential with strong coupling asymptotics 

y(A)~A^ as A^oo, (6.18) 
taking into account ( |6.13 ), equation ( |6.17 ) is modified to 



c 



s 



small 327r 



as rh ro. (6.19) 



where rg is the position of the singularity in the zero temperature solution. 

For confining theories, the limit r^ — )• ro corresponds to T — )• oo on the small BH 
branch. However, one can show that the result ( |6.19| ) holds quite generally, regardless of 
whether the zero T theory confines or not. In particular, for the non-confining theories — 
that is either when Q < 4/3 or when Q = 4/3 but the subleading term in the potential 
vanishes at the singularity — there is only the large black-hole branch and the limit r^ — >• r^ 
corresponds to the zero T limit of this BH. Thus, we also learn that there exist holographic 
models that correspond to non-confining gauge theories whose zero T limit yield a constant 
C/s. This constant approaches zero as Q — )• 0, i.e. in the limiting AdS case. 

We also see that the asymptotic value of C/s in the small BH branch is close to the 
value of C/s near Tc- We shall give an explanation of this fact in the next subsection. Using 
the asymptotic formula (6.19), the fact that Q > | for confinement and Q < for the 



IR singularity to be good and repulsive we may obtain a range of values where we expect 
C/s to vary, namely 



± < i 

6tt s 



< i-. (6.20) 

small, asymptotic 



A final observation concerns the coefficient ci,{Xh) in ( |6.13| ). This part is the only 
input from the solution of the fiuctuation equation, the rest of ( |6.13| ) is fixed by the dilaton 
potential entirely. We plot the numerical result for in fig |l^ as a function of the coupling 
at the horizon A/^. 

First of all, Figure 15 provides a check that, the approximate bound of [^] \ch\ > 1, 



is satisfied in the entire range. One also observes Cf, approaches to 1 in the IR and UV 
asymptotics. These facts can be understood analytically: In the UV (near the boundary) 
it is because of the boundary condition Cb = 1. In the IR, it is more subtle, and it is 
explained in appendix B of |f6|. 



^See the discussion at appendix B of fF^ . 
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Figure 15: The coefScient \cb\ of equation (3.13) as a function of Xh 



Finally, we observe that the deviation of C5 from the asymptotic value 1 is maximum 
around the phase transition point Ac. In fact, we numerically observed that the top of 



the curve in figure 15 coincides with Ac to a very high accuracy. Whether this is just a 



coincidence or not, remains to be clarified. 
6.3 The adiabatic approximation 

Motivated by the Chamblin-Reall solutions |77], Gubser et al. [|l6| proposed an approximate 
adiabatic formula for the speed of sound. In the case when V' /V is a slowly varying function 
of (/?, proposes the following formulae for the entropy density and the temperature: 

8 /■'^h Y 

logs = --J + ..., (6.21) 

/ 1 y 8v 



io,T = J (6.22) 

where the ellipsis denote contributions slowly varying in tph'^'^- 

It is very useful to reformulate this approximation using the method of scalar variables, 
which in turn allows us to extract the general T dependence of most of the thermodynamic 
observables in an approximate form. Here, we apply this formalism to the computation of 
C/s. The method of scalar variables and the details of the adiabatic approximation in the 



scalar variables are given in |76|. 



For the scalar variable X = ^7 the adiabatic approximation means 

(0.23, 



The fluctuation equation ( |6.7| ) greatly simplifies with ( |6.23 ). In fact, as shown in |76], the 



solution becomes independent of (p. With unit normalization on boundary, the adiabatic 



^^Various coefficients in these equations differ from due to our different normalization of the dilaton 
kinetic term. 
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Figure 16: Comparison of the exact C/s with the adiabatic approximation in the variable Xh- 
Sohd(red) curve is the full numerical result and the dashed(blue) curve follows from ( 6.24 ). 




T/Tc 



Figure 17: Comparison of the exact (/s with the adiabatic approximation in variable T. Solid(blue) 



curve is the full numerical result and the dashed(red) curve follows from (6.24) 



solution in the entire range of 9? £ {— 00, iph) becomes hadbi^) 
coefficient Cf, in ( |6.13 ) becomes unity, hence: 



c 



3 ( V'{^h) \ 
327r V V{^h) J 



1. Consequently, the 



.24) 



We plot this function in in figure |T^, where we also provide the exact numerical 



result for comparison. Note that in figure |16| the whole large black-hole branch has been 
compressed at the left of the figure for Xh ^ 0.04 The same functions in the variable T/Tc 



are plotted in figure |17. 

The validity of the adiabatic approximation equation ( |6.23 ), is determined by the rate 
at which V /V varies with In particular, the approximation becomes exact in the limits 
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where V' jV becomes constant. This happens for a constant potential or a potential that 
is a single power of A (exponential in This is the case in the UV (99 — t- —00, where the 
potential becomes a constant) and the IR (ip — >• +00 where the potential becomes a power 
law.) . Therefore equation ( |6.24 ) allows us to extract the analytic behavior of Qjs in the 



limits — )■ ibcxD. 

The numerical values one obtains from ( |6.24 ) in the intermediate region may differ 



from the exact result ( 6.13| ) considerably, especially near Tc- However, we expect that the 



general shape will be similar. 

Finally, the adiabatic approximation hints at why, in the particular background that 
we study, C/s at Tc is close to the limit value ( |6.17 ): In order to see this we rewrite ( 6.24| ) 
as 

^ = (6.25) 

In the limit ( |6.17 ) we have X — )■ —1/2. The only other point where X = —1/2, is at 
the minimum of the string frame scale factor (p^. This is the point where the confining 
string saturates [^4|. On the other hand, we expect on general physical grounds that the 
de-confinement phase transition happens near this point, i.e. (/Jc ~ V'*- Thus, the adiabatic 
formula predicts that C/s{ipc) be close to the limit value l/Gvr.^^ 

6.4 Buchel's bound 



In [78 1, Buchel proposed a bound for the ratio of the bulk and shear viscosities, motivated 
by certain well- understood holographic examples. In 4 space-time dimensions the Buchel 
bound reads, 

i>2(i-cj). (6.26) 

We note that the bound is proposed to hold in the entire range of temperature from Tc to 
00. This bound is trivially satisfied for exact conformal theories such as A/" = 4 YM, and 



saturated in theories on Dp branes |7g, With the numerical evaluation at hand, we 
can check (|]26|) in our case. In figure ^ (a) we plot the LHS and RHS of the bound 
We clearly see that the bound is satisfied for all temperatures. As expected, both the LHS 
and the RHS of ( |6.26| ) vanishes in the high T conformal limit. 

A clear picture of Buchel's bound is obtained by defining the function: 

in terms of which the bound is simply C > 1. In Figure (b) we show the function 
C{T) obtained numerically in our IHQCD model, between Tc and 5Tc. The values of this 
function are mildly dependent on temperature, and are between 1.5 and 2, the same range 
of values that were recently considered in the hydrodynamic codes by Heinz and Song pO|. 



^^This argument may break down for two (dependent) reasons: First of all the adiabatic approximation 
becomes lees good near ipc- This is because, in this region V' /V varies relatively more rapidly as a function 
of tfi. Secondly, precisely because of this, even though ipc is not far away from ip, the difference can result 



in a considerable change in the value of ("/s through (3.24). 

^^Since this theory contains two derivatives only, ^ has the universal value l/47r. 
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Figure 18: (a) Comparison of C^jr] (solid line) and the RHS of (6.26) (dashed line), obtained using 
the speed of sound of the IHQCD model j22|. (b) Plot of the function C(r) defined in equation 
( |6.27 ) as a function of temperature. The horizontal dashed line indicates where Buchel's bound is 
saturated. We see that the bound is satisfied in the entire range of temperatures. 



We may also investigate the fate of the bound at large T. using the asymptotics of C/s 
in (l6l5|) 

C 1 1 

+ • • • , as T ^ oo. 



27^1og^(^r/A) 



and 



1 _ 4 1 326logil°giir 
3 — — 7 ^ h 



9 



log^ 



+ 



from [M, that can be rewritten as 



.28) 



.29) 



3 ' 



1 32&log(log(7l_ 

+ — ^ — r X + 



log^ 



Tc 



(6.30) 



where b = 



3-34 
2-121 



is the ratio of the two-loop to the one-loop squared /3-function 



coefficients in large- A'^c YM. 

Since in this class of models rj/s = l/Air exactly we obtain 



lim , , 

T^oo 2(1/3- c2) 



(6.31) 
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in agreement with a recently derived general formula, in Einstein dilaton gravity, |]81| , 

lim = 2vri^ cot ( ^) (6.32) 

T^oo 2(1/3 - c2) 4-2A V 4 y 

where A is the scaling dimensions of the scalar operator in the UV, that is marginal in our 
case. 



It has also been suggested recently, [p2, 83|, that the speed of sound squared, in Einstein 



dilaton gravity asymptotes to 1/3 at high temperatures from below. This is evident in our 
asymptotic formula ( |6.30 ), although the formulae in |82, 83] fail to capture correctly the 



marginal case that is relevant here. 

7. The drag force on strings and heavy quarks 

We will now consider an (external) heavy quark moving through an infinite volume of gluon 
plasma with a fixed velocity v at a finite temperature T The quark feels a drag 

force coming from its interaction with the plasma and an external force has to be applied 
in order for it to keep a constant velocity. In a more realistic set up one would like to 
describe the deceleration caused by the drag. 

The heavy external quark can be described by a string whose endpoint is at the bound- 
ary. One can accommodate flavor by introducing D-branes, but we will not do this here. 
A first step is to describe the classical string "trailing" the quark. 

We consider the Nambu-Goto action on the world-sheet of the string. 

Sng = J dadr^det {-gMNdaX^ dpX^) , (7.1) 

where the metric is the string frame metric. The ansatz we are going to use to describe 



the trailing string is, |11] 



X^=vt + i{r), X2=X3 = 0, (7.2) 

along with the gauge choice 

cr = r, T = t , (7.3) 

where r is the (radial) holographic coordinate. The string is moving in the X^ direction. 

This is a "steady-state" description of the moving quark as acceleration and decel- 
eration are not taken into account. For a generic background the action of the string 
becomes ^ 

S = J ^^^^ V-9oo9rr - googiiS,'^ - gnOrrV^ ■ (7.4) 

Note that goo is negative, and we should check whether our solution produces a real action. 
For example a straight string stretching from the quark to the horizon is a solution to the 
equations of motion but has imaginary action. 

We note that the action does not depend on ^ but only its derivative, therefore the 
corresponding "momentum" is conserved 

^5 = ~ 27r/2 / 2 • (^-^^ 
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We solve for ^' to obtain 



■googrr - gngrrV 



Vgoogii (1 + 5'oo5'ii/(2vr£27rg)2 



(7.6) 



The numerator changes sign at some finite value of the fifth coordinate r <j . For the solution 
to be real, the denominator has to change sign at the same point. We therefore determine 
rg via the equation 

goo{rs) + gii{rs)v'^ = , (7.7) 



and the constant momentum 



Writing the string-frame metric as 



^^7^) becomes 



gooirs)gn{rs 
(2vr£2)2 



— / dt^ + dx ■ dx 



= firs) 



(7.8) 

(7.9) 
(7.10) 



The induced world-sheet metric is therefore 



p2A.(r) 



f(r) Y e4^s('-)/(r)-e4'4s(rs)„2 

4As(r) j2(-j.^_^4g4As(rs) 



7(7) y e^^s(r)f(r)-e't^s(rs)v^ p (r) [e^As (r) /(r) - t,2 e4^s (rs ) J 



(7.11) 

We can change the time coordinate to obtain a diagonal induced metric t = t + ((r) 



with 



,2A4r,) 2 



f{r)^/U{r) - v2)(e4>l.W/(r) - e4^4r.)^;2) 



The new metric is 

ds^ = e^^^^'-) 



-(/(r)-t;^)dr^ + 



dr^ 



(7.12) 



and near r = it has the expansion 



ds^ 



-f{r,)e^^^^''^\r-r,) + 0{{r-rsf) 



dT^+ 



,2A4rs) 



{4.v-'A>,{r,)+ f'{rMr-rs) 



+ 0{l) 
(7.13) 



dr^ 



This is a world-sheet black-hole metric with horizon at the turning point r = rg. 
7.1 The drag force 

The drag force on the quark can be determined by calculating the momentum that is lost 
by flowing along the string into the horizon: 



drag 



dpi 



1 500511?' 



dt 27r£2 



(7.14) 
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This can be obtained by considering the world-sheet Noether currents and expressing 
the loss of momentum as APJ'^ = J" H^. This may be evaluated at any value of r, but it is 
more convenient to evaluate it at r = r^. 
We finally find that 



drag 



~^V-9ooirs)gii{rs) ■ 



(7.15) 



Using the form ( |7.9D of our finite-temperature metric in the string frame we finally 
obtain 



rag 



27r£2 



27r£2 



(7.16) 



where in the second equality we expressed the force in terms of the Einstein-frame scale 
factor and the "running" dilaton. Substituting from ( [7.10 ) we obtain 



rag 



27r£2 



(7.17) 



Before proceeding further, we will evaluate the drag force for the conformal case of = 4 
SYM where 



f{rs) = 1 - (vrTrs 



Substituting in ( [7.17| ) we obtain, @-|l3, 

i^conf = r2 



(7.18) 



(7.19) 



Moving on to YM, to compute the drag force from equation (7.17) we must first 
determine is in the IHQCD model. In this setup there is no analog of the = 4 SYM 
relation ( |7.18| ) between £, ig and A. Rather, the fundamental string length ig is determined 
in a bottom-up fashion, by matching the effective string tension to the QCD string tension 
CTc derived from the lattice calculations. We obtain 



(7.20) 



where r^, is the point where the zero-temperature string scale factor (at T=0) Ag^oif) has 
a minimum. For a typical value of ac ~ (440 MeV)"^ ||6l[] we find 



L = 6.4 



(7.21) 



where £ is the radius of the asymptotic AdS space. 

On the other hand, unlike in = 4 SYM, in the IHQCD model the value of the 



coupling X{rs) in equation ( 7.17 ) is not an extra parameter to be fixed by hand, but rather 
it is determined dynamically together with the background metric. 
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Figure 19: In this figure tlie ratio of the drag force in improved holographic QCD to the drag force 
in M — 4: SYM is shown. The ratio is computed for different velocities as a function of temperature. 
The 't Hooft coupling for the JV = 4 SYM theory is taken to be 5.5. We chose this value as it is 
considered in the central region of possible values for the 't Hooft coupling. It is seen that as the 
velocity increases the value of the ratio decreases. 
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Figure 20: In this figure the ratio of the drag force in improved holographic QCD to the drag force 
in A/" = 4 SYM is shown. The ratio is computed for different temperatures as a function of velocity. 
The 't Hooft coupling for the Af = 4 SYM theory is taken to be 5.5. As temperature increases the 
value of the ratio decreases. 



7.2 The relativistic asymptotics 

When f — )■ 1 then — )• and we approach the boundary. Near the boundary (r — t- 0) we 
have the following asymptotics of the scale factor and the 't Hooft coupling, 111] 



f{r) ^ 1- 
and 



l + O 



A(r) 



log(Ar) 



Air) 



1 



Po log(rA) 
where r/j is the position of the horizon. 

We therefore obtain for the turning point 



+ 0(log(rA)-2) 



l + O 



1 



log(Ar) 



^3(1 



l + O 



1 



log(l — v'^ 



A(r.) 



(3o log [1 - v"^' 



+ 



+••• 

(7.22) 

(7.23) 



(7.24) 
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and the drag force 

\/TTaHrH)\i{r,) v 



We also use 



47rM3 Ar2 Ar2 



where s(T) the entropy per unit three- volume, and we write the relativistic asymptotics of 
the drag force as, 



^drag ^ -^^^^ ;T + • • • (7-27) 



/45 Ts{T) V 



[-^\og [1-^2] 



£2 y 4iv2 

The force is proportional to the relativistic momentum combination v / \/\ — modulo 
a power of log [l — v^] . This factor is present because, as argued in the asymptotic 
metric is AdS in the Einstein frame instead of the string frame. Its effects are not important 
phenomenologically. 

7.3 The non-relativistic asymptotics 

We now consider the opposite limit, f — t- 0. In this case the turning point asymptotes to 
the horizon, — >• and we have the expansion 



f{T)^\T:T{rh-r) + 0{(Th-rf) , r, = r;, - £^ + (7.^ 



2 



and 



2-kT ' SttT \{rh) 



(7.29) 



where primes are derivatives with respect to the conformal coordinate r. 
7.4 The diffusion time 

For a heavy quark with mass Mq we may rewrite ( [7. 191) as 

"conf = — = -zP ' P= 2 ^ ' 



dt ' ^ 



where the first equation defines the diffusion time r. In the conformal case, the diffusion 
time is constant, 

Tconf = ^ 



vr 



T2 
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This is not anymore the case in QCD, where r defined as above is momentum dependent. 
We may still define it as in ( [7.30| ) in which case we obtain the following limits 



4iV2 



p—^oo 



2X1 



+ 



(7.32) 







lim 

p — >\j 
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Figure 21: In this figure the ratio of the diffusion time in the Improved Holographic QCD model 
to the diffusion time in A/" = 4 SYM is shown. The 't Hooft coupling for A/" = 4 SYM is taken to 
be A = 5.5. The heavy quark has a mass of Mq = 1.3GeV. Note that with the definition of the 



diffusion time in (7.30) the ratio is the inverse of the ratio of the forces. A similar plot is valid 
for the bottom quark as well, as the mass drops out of the ratio, although the energy scales are 
different. In this plot the x-axis is taken to be in MeV units. As temperature increases the ratio 
also increases. 



7.5 Including the medium-induced correction to the quark mass 

In order to estimate the diffusion time of a quark of finite rest mass, we must take into 
account the fact that the mass of the quark receives medium-induced corrections. In 
other words, the mass appearing in equation ( 7.30| ) is a temperature-dependent quantity, 
Mq{T) 7^ Mq{T = 0). The ratio Mq{T)/Mq{0) can be estimated holographically by repre- 
senting a static quark of finite mass by a static, straight string 2'' stretched along the radial 
direction starting at a point r = Vq ^ 0. At zero temperature, the IR endpoint of the string 
can be taken as the "confinement" radius, r^,, where the string frame metric reaches its 
minimum value; At finite temperature, the string ends in the IR at the BH horizon^^. The 



^^This representation ignores the fact that the kinetic mass of a moving quark may be different from the 
static mass 

^*It would stop at the confinement radius if the latter were closer to the boundary than the horizon, 
i.e. if r,{T) < rniT). However, in the model we are considering, in the large BIf branch we find that the 



relation th < r* is always satisfied. 
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masses of the quark at zero and finite T are related to the world-sheet action evaluated on 
the static solution {t = t,a = r) : 



M,(0) 



rh 



dre2^MA4/3(r). (7.34) 



The value can be fixed numerically by matching Mq{0) to the physical quark mass, 
and translating the fundamental string tension in physical units by using the relation (|7.2C| ), 
with cjc = {AAOMeVy. This makes Mg{T) a function of Mg{0). The ratios Mg{T)/Mg{0) 



we found numerically in the model under consideration is shown in figure 22 for the Charm 
(M(0) = l.bGeV) and Bottom (M(0) = i.bGeV) quarks. The fact that, in the deconfined 
plasma, the quark mass decreases with increasing temperature is a direct consequence of 
the holographic framework^^, since for higher temperature, the distance to the horizon 
is smaller. An indication that this result may be in the right direction comes from the 
lattice computation of the shift in the position of the quarkonium resonance peak at finite 



temperature |85|: in the deconfined phase the charmonium peak moves to lower mass at 
higher temperature. Our result for the medium-induced shift in the constituent quark mass 
is consistent with these observations. 
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Figure 22: Ratios between the thermal mass and the rest mass of the Charm (curve labeled "c") 
and Bottom (curve labeled "b" ) quarks, as a function of temperature. 



We can now write the diffusion time from eqs. (7.17) and ( 7.3C| ) as: 

Mg{T) /Ao(n)^^/3 



t{T,v) 



,2Ao{r^)-2A{rs) 



(7.35) 



where once again we have eliminated the fundamental string length using equation ( 7.2C| ) . 
Given a set of zero- and finite-temperature solutions, equation ( [7.35| ) can be evaluated 



'For a possible field theoretical explanation of this phenomenon, see |84 
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numerically for different values of the velocity and different quark masses. The results for 
the Charm (Mg(0) = 1.5 GeV) and Bottom (M = 4.5 GeV) quarks are displayed in figure 
23|. 
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Figure 23: Diffusion time for the Charm and Bottom quarks, as a function of energy, for different 
ratios of the temperature to the IHQCD transition temperature Tc- 

7.6 Temperature matching and diffusion time estimates 

An important question is how we should choose the temperature in our holographic model 
in order to compare our results with heavy-ion collision experiments. This is nontrivial, 
since our setup is designed to describe pure SU{Nc) Yang- Mills, whereas at RHIC temper- 
atures there are 3 light quark flavors that become relevant. As a consequence, the critical 
temperatures and the number of degrees of freedom of the two theories are not the same: 
for pure SU (iVc) Yang Mills we have — 1 degrees of freedom and a critical temperature 
around 260 MeV; For SU (A^c) QCD with Nf flavors the number of degrees of freedom is 

— 1 -|- NcNf, and the transition temperature is lower, around 180 MeV. 

In IHQCD, the transition temperature in physical units was calculated to be Tc = 
247 MeV |^2[, i.e. close to the lattice result for the pure YM deconfining temperature. 
From now on, this is the value we will mean when we refer to Tc- This is also close to 
the temperature of QGP at RHIC, which we will denote Tqgp, and is estimated to be 
around 250 MeV. Since this value is uncertain, below we give our results for a range of 
temperatures between 200 MeV and 400 MeV. The higher temperatures will be relevant 
for the LHC heavy-ion collision experiments (see e.g. [^]). 

Based on these considerations, there are different ways of fixing the temperature (see 
e.g. the recent review |3^): one direct and two alternative schemes (that we call the energy 
and entropy scheme). 

• Direct scheme: The temperature of the holographic model is identified with the 
temperature of the QGP in the experimental situation (at RHIC or LHC), ~ 




Energy scheme: One matches the energy densities, rather than the temperatures. 
The energy density at RHIC is approximately (treating the QCD plasma as a free 



gas^o.) eQGP ^ (7rVl5)(iV2 _ ^ ^ NcNf){TQGp)^. For N^ = Nf = 3, asking that our 
energy density matches this value requires us to consider the holographic model at 
temperature T^^^^^ given by 

^^Hcd{Ti^^ ^ ll.2{TQGpf (7.36) 



Entropy scheme: Instead of matching the energy densities, alternatively one can 
match the entropy density s, which for the QGP, in the free gas approximation, is 
given by (Tqgp ^ 47r^/A5{N^ - 1 + NcNf){TQGp)^. This leads to the identification: 

^ihUTlnlj = ^^-HTQGpf (7.37) 



Tqgp (MeV) 


Tqgp/Tc 


^Scd (MeV) 


rp{e) Jrp 

^ihqcd/ 


T^ca (MeV) 


^ihqcdl 


190 


0.77 


259 


1.05 


274 


1.11 


220 


0.89 


290 


1.18 


302 


1.23 


250 


1.01 


325 


1.31 


335 


1.35 


280 


1.13 


361 


1.46 


368 


1.49 


310 


1.26 


398 


1.61 


402 


1.63 


340 


1.38 


434 


1.76 


437 


1.77 


370 


1.50 


471 


1.90 


472 


1.91 


400 


1.62 


508 


2.06 


507 


2.05 



Table 3: Translation table between different temperature identification schemes. The first two 
columns display temperatures in the direct scheme, (in which the temperature of the holographic 
model matches the physical QGP temperature) and the corresponding ratio to the IHQCD critical 
temperature, that was fixed by YM lattice results at = 247 MeV The third and fourth 

columns display the corresponding temperatures (and respective ratios to Tc) in the energy scheme, 
and the last two in the entropy scheme. 

The temperature translation table between the various schemes is shown in Table |3|. 
In that table, = 247MeV is the deconfining temperature of the holographic model. 

In Figure ^ we show the comparison between the diffusion times, as a function of 
initial quark momentum, in the different schemes for the Charm and Bottom quarks, at 
the temperature Tqgp = 250MeV. 

The results for the diffusion times at different temperatures, computed at a reference 
heavy quark initial momentum p « 10 GeV , are displayed in Tables |^ and ^. We see that 
there is little practical difference between the entropy and energy schemes; on the other 
hand the difference between the direct scheme and the two alternative schemes can be quite 
substantial. 

^"Tfiis is itself an approximation, since as we know both from experiment and in our holograpliic model, 
the plasma is strongly coupled up to temperatures of a few Tc 
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Figure 24: Diffusion times for the Charm and Bottom quarks, as a function of initial momentum, 
at Tqgp — 250 MeV. The different lines represent the in the direct scheme (solid), energy scheme 
(dashed) and entropy scheme (dash-dotted) , all corresponding to the same temperature Tqq p = 250 
MeV. 



Tqgp, MeV 


Tcharm (fm/c) 

(direct) 


Tcharm (fm/c) 

(energy) 


Tcharm (fm/c ) 

(entropy) 


220 




4.0 


3.6 


250 


5.7 


3.1 


3.0 


280 


4.3 


2.6 


2.5 


310 


3.5 


2.1 


2.1 


340 


2.9 


1.8 


1.8 


370 


2.5 


1.5 


1.5 


400 


2.1 


1.3 


1.3 



Table 4: The diffusion times for the charm quark are shown for different temperatures, in the 
three different schemes. Diffusion times have been evaluated with a quark initial momentum fixed 
at p w 10 GeV. 



8. Jet quenching parameter 

In this Section we discuss the jet quenching parameter in the class of holographic models 
under consideration, and we estimate its numerical value for the concrete model with 
potential ( [4.1D and parameters fixed as in [^2|. For the holographic computation, we will 
follow p9|, |40|]. There is another method available |45], but we will not use it here. 



The jet-quenching parameter q provides a measure of the dissipation of the plasma 
and it has been associated to the behavior of a Wilson loop joining two light-like lines. We 
consider two light-like lines which extend for a distance L~ and are situated distance L 
apart in a transverse coordinate. Then q is given by the large L"*" behavior of the Wilson 
loop 

W^e~^^^~^\ (8.1) 
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10.1 
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6.1 
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8.6 
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5.3 


340 


7.5 


4.7 


4.7 


370 


6.6 


4.1 


4.1 


400 


5.8 


3.6 


3.6 



Table 5: Diffusion times for the bottom quark are shown for different temperatures, in the three 
different schemes. Diffusion times have been evaluated with a quark initial momentum fixed at 
p w 10 GeV. 



We consider the bulk string frame metric 



(8.2) 



To address the problem of the Wilson loop we make a change of coordinates to light cone 
coordinates for the boundary theory 



X = Xl + t X = Xl —t 



(8.3) 



for which the metric becomes 



ds^ = e^^^ (^dxl + dxl + ^{l- f){dx\ + dx'i) + {l + f)dx+dx- + . (8.4) 

The Wilson loop in question stretches across X2-, and lies at a constant X-\.,x^. It is conve- 
nient to choose a world-sheet gauge in which 



= r, X2 = cr . 

Then the action of the string stretching between the two lines is given by 

S = ^J dadT^-det{gMNdaX^dpX^) 
and assuming a profile of r = r{a) we obtain 

The integrand does not depend explicitly on X2, so there is a conserved quantity, c: 

'—-r = — 
dT> V2 



(8.5) 



(8.6) 



(8.7) 



(8.8) 
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Figure 25: In this figure the combination (1 — f)e^^'' is plotted as a function of the radial distance, 
for several temperatures. The radial distance is given in units of the horizon position rc for the 
black hole at the critical temperature Tc. All curves stop at the corresponding horizon position. 



which leads to 

'e4^^(l-/) 



r" = f[ \ - 1 • (8-9) 



A first assessment of this relation involves determining the zeros and the region of positivity 
of the right-hand side. / is alvi^ays positive and vanishes at the horizon. For the second 
factor we need the asymptotics of e^^''{l — /). This factor remains positive and bounded 
from below in the interior and up to the horizon. It vanishes however logarithmically near 
the boundary as 



8 

This is unlike the conformal case where we obtain a constant 



1 



log(Ar) / 



(8.10) 



, =i^Tir (8.11) 

coniormal 

Because of this, for fixed c, there is a region near the boundary where r'^ becomes 
negative. At this stage we will avoid this region, by using a modified boundary at r = e. 
We will later show that this gymnastics will be irrelevant for the computation of the jet 
quenching parameter, as it involves effectively the limit c — t- 0. 

We will place the modified boundary r = e a bit inward from the place r = rmin where 
the factor \ ■" — 1 vanishes: 

e^^'^'-'-Hl-f{rmin)) = c'' (8.12) 

Therefore we choose Tmin < £• 
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Then, in the range e < r < the factor - — — 1 is positive for sufficiently small 
c. In this same range, r' vanishes only at r = rh- This is the true turning point of the 
string world-sheet. This is also what happens in the conformal case. It is also intuitively 
obvious that the relevant Wilson loop must sample also the region near the horizon. 

The constant c is determined by the fact that the two light-like Wilson loops are a 
X2 = L distance apart. 

- = / = . (8.13) 

The denominator vanishes at the turning point. The singularity is integrable^^ . Therefore, 
as we are interested in the small L region, it is obvious from the expression above that that 
c must also be small in the same limit. 

This relation can then be expanded in powers of c as 

2c A y7(W) 2 V/(l-/)3 
Therefore to leading order in L 

+ 0{L^) (8.15) 



2/. 



We are now ready to evaluate the Nambu-Goto action of the extremal configuration 
we have found. Starting from ( |8.7| ), we substitute r' from ( |8.9| ), and change integration 
variable from X2 — )• r to obtain 

- '^-r - , ^y-'^ (8.16) 



As in p9|, 40 1, we subtract from equation (p. 16) the action of two free string straight 



world-sheets that hang down to the horizon. To compute this action a convenient choice 
of gauge is x_ = r, r = a. The action of each sheet is 



The subtracted action is therefore: 



5 _ 2^0 = ^ / , + O(c^) , (8.18) 



Using now ( p.l5| ) to substitute c we finally obtain 



Even if we choose e — rmin, the new singularity a± r — r,nin is also integrable as suggested from (8.10) 
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So far we have evaluated the relevant Wilson loop in the fundamental representation 
(by using probe quarks). On the other hand, the Wilson loop that defines the jet-quenching 
parameter is an adjoint one. We can obtain it in the large- A^c limit from the fundamental 
using tr Adjoint = ^^Fundamental- finally extract the jet-quenching parameter as 



1: 



We are now ready to remove the cutoff. As the integral appearing is now well-defined up 
to the real boundary r = we may rewrite it as 

/ e ^^^^ ^ ^^^^ ^ I e^^^ ar ^^^^^^ 

In [|76| we obtain the small e estimate of 1(e) that vanishes as ~ e(loge)3. We may 
finally write'^'^ 

. \/2 1 



fh dr 

Jo 



.22) 




Figure 26: In this figure the ratio of the jet quenching parameter in our model to the jet quenching 



parameter in A/" = 4 is shown. The integral present in equation (8.2C) has been numerically 
calculated from an effective cutoff at r = r^i/lOOO. The jet quenching parameter in A/" = 4 SYM 
has been calculated with X'tHooft = 5.5. 



From equation ( |8.22| ) we obtain, in the conformal case: 



r [-1 

gconformal = M V2X (8.23) 

The conformal value, for the median value of A = 5.5 and T ~ 250 MeV gives (?conformai — 
1.95 GeV^/fm where we used the conversion 1 GeV~ 5 fm~^. 



^^In practise, the previous discussion including regularizing the UV is academic. The numerical calcula- 
tion is done with a finite cutofT where the boundary conditions for the couplings are imposed. 



-52- 



q, GeV /fm 
40 - 

30 

20 

10 



/ 


ihQCD 






// 






ihQCD bottom 








ihQCD charm 




AdS, 


T/Tc 




2.5 3.0 





Figure 27: The jet quenching parameter q for the Improved Holographic QCD model and J\f = 4 
SYM is shown in units of GeV^/ fm for a region close to T = Tc- The smallest dashed curve is the 
ihQCD result with an effective cutoff of r cutoff ~ rji/1000. The small dashed curve is the ihQCD 
result with the cutoff from the mass of the Bottom quark. The medium dashed curve has a cutoff 
coming from the Charm mass and and largest dashed curve is the Af — A SYM result. 






ihQCD 




AdSj 



T/Tc 



Figure 28: The jet quenching parameter q for the Improved Holographic QCD model (lower curve) 
and Af — 4: SYM (upper curve) are shown in units of GeV^ / fm for temperatures up to T = ATc. 
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Table 6: This table shows the jet quenching parameter q computed with different cutoffs for the 
different temperatures shown in the first column. The computation is done in the direct scheme. 
The second column shows q with a cutoff at r cutoff — ''h/lOOO, where is the location of the 
horizon. In accordance with the conclusions of appendix q does not change significantly as we vary 
the cutoflF from r,j/1000 to r^/lOO. 



Numerical evaluation of equation ( ^.22| ) in the non-conformal IHQCD setup^^ gives 
us a value of q which is lower (at a given temperature) than the conformal value, as 



In this case, the value of Is appearing in equation (8.22) is fixed as explained in Section 4 



-53- 



l^^^r, A/fcT/ 

J^QGP, iviev 


q [LjeV /jm) 


q [Lrev / jm) 


9 i<-rev/ / jm) 




I U-il cUL J 


energy J 


(^enriopyj 


990 




n Q 


1 n 


9^n 




1 9 


1 

± .o 


280 


0.8 


1.6 


1.7 


310 


1.1 


2.1 


2.2 


340 


1.4 


2.7 


2.8 


370 


1.8 


3.4 


3.4 


400 


2.2 


4.2 


4.2 



Table 7: This table displays the jet quenching parameter q using the three different comparison 
schemes. For lower temperatures the "entropy scheme" gives higher values. As energy is increased 
the energy and entropy schemes temperatures start to coincide and there is little difference in the 
jet quenching parameter as well. 
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Table 8: This table displays the jet quenching parameter q using the three different comparison 
schemes with an effective cutoff provided by the mass of the Charm quark. Again, for lower 
temperatures the "entropy scheme" gives higher values. As energy is increased the energy and 
entropy schemes temperatures start to coincide and there is little difference in the jet quenching 
parameter as well. Also when the temperature approaches the quark mass the picture of the heavy 
quark as a hanging string collapses and results are not reliable. 



shown in Figures 27 and 28. Tables ^ to ^ display the numerical values of the jet 
quenching parameter at different temperatures in the experimentally relevant range, in 
different temperature matching schemes. 



9. Discussion and outlook 



The construction presented in this paper offers a holographic description of large- A'^c Yang- 
Mills theory that is both realistic and calculable, and in quite good agreement with a large 
number of lattice results both at zero and finite temperature. 

It is a phenomenological model; as such it is not directly associated to an explicit string 
theory construction. In this respect it is in the same class as the models based on pure 
AdS backgrounds (with hard or soft walls) 0, |^ or on IR deformations of the AdS 

metric [8£, 90]. In comparison to them however, the IHQCD approach has the advantage 
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Table 9: This table displays the jet quenching parameter q using the three different comparison 
schemes with an effective cutoff provided by the mass of the Bottom quark. The results are close 
to the q results computed in Table |^ since the mass of the Bottom quark is much larger than the 
temperatures we examine. 

that the dynamics responsible for strong coupling phenomena (such as confinement and 
phase transitions) is made explicit in the bulk description, and it is tied to the fact that the 
coupling constant depends on the energy scale and becomes large in the IR. This makes 
the model consistent and calculable, once the 5D effective action is specified: the dynamics 
can be entirely derived from the bulk Einstein's equation. The emergence of an IR mass 
scale and the finite temperature phase structure are built-in: they need not be imposed 
by hand and do not suffer from ambiguities related to IR boundary conditions (as in hard 
wall models) or from inconsistencies in the laws of thermodynamics (as in non-dynamical 
soft wall models based on a fixed dilaton profile |^] or on a fixed metric p9| , pO[]). More 
specifically, in this approach it is guaranteed that the Bekenstein-Hawking temperature 
of the black hole matches the entropy computed as the derivative of the free energy with 
respect to temperature. 

With an appropriate choice of the potential, a realistic and quantitatively accurate de- 
scription of essentially all the static properties (spectrum and equilibrium thermodynamics) 
of the dynamics of pure Yang-Mills can be provided. The main ingredient responsible for 
the dynamics (the dilaton potential) is fixed through comparison with both perturbative 
QCD and lattice results. It is worth stressing that such a matching on the quantitative 
level was only possible because the class of holographic models we discuss generically pro- 
vides a qualitatively accurate description of the strong Yang-Mills dynamics. This is a 
highly non-trivial fact, that strongly indicates that a realistic holographic description of 
real-world QCD might be ultimately possible. 

Although the asymptotics of our potential is dictated by general principles, we base our 
choice of parameters by comparing with the lattice results for the thermodynamics. There 
are other physical parameters in the 5D description that do not appear in the potential: 
the 5D Planck scale, that was fixed by matching the free field thermodynamics in the limit 
T — 7- oo; the coefficients in the axion kinetic term, that were set by matching the axial 
glueball spectrum and the topological susceptibility (from the lattice). The quantities that 
we use as input in our fit, as well as the corresponding parameters in the 5D model, are 
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shown in the upper half of Table |[ 

The fact that our potential has effectively two free parameters depends on our choice 
of the functional form. This functional form contains some degree of arbitrariness, in that 
only the UV and IR asymptotics of V{X) are fixed by general considerations (matching 
the perturbative /^-function in the UV, and a discrete linear glueball spectrum for the 
IR). Therefore the results presented in this paper offer more a description, rather than a 
prediction of the thermodynamics. 

Nevertheless, there are several quantities that we successfully "postdict" (i.e. they 
agree with the lattice results) once the potential is fixed: apart from the good agreement of 
the thermodynamic functions over the whole range of temperature explored by the lattice 
studies (see Figures |5|, and P), they are the lowest glueball mass ratios and the value 
of the critical temperature. The comparison of these quantities with the lattice results 
is shown in the lower half of Table ^ and one can see that the agreement is overall very 
good. Moreover the model predicts the masses of the full towers of glueball states in the 
0+~,0++,2++ families. 

The fact that IHQCD is consistent with a large number of lattice results is clearly not 
the end of the story: its added value, and one of the main reasons for its interest lies in its 
immediate applicability beyond equilibrium thermodynamics, i.e. in the dynamic regimes 
tested in heavy-ion collision experiments. This is a generic feature of the holographic 
approach, in which there are no obstructions (as opposed to the lattice) to perform real- 
time computations and to calculate hydrodynamics and transport coefficients. IHQCD 
provides a framework to compute these quantities in a case where the static properties 
agree with the real-word QCD at the quantitative level. Therefore the bulk viscosity, drag 
force and jet quenching parameters were computed in IHQCD: 

Bulk Viscosity: The bulk viscosity was computed by calculating the low frequency 
asymptotics of the appropriate stress tensor correlator holographically. The result is that 
the bulk viscosity rises near the phase transition but stays always below the shear viscosity. 
It fioats somewhat above the Buchel bound, with a coefficient of proportionality varying 
between 1 and 2. Therefore it is expected to affect the elliptic fiow at the small percentage 



level [g^, 33]. Knowledge of the bulk viscosity is important in extracting the shear viscosity 
from the data. This result is not in agreement with the lattice result near Tc. In particular 
the lattice result gives a value for the viscosity that is ten times larger. 

The bulk viscosity keeps increasing in the black-hole branch below the transition point 
until the large BH turns into the small BH at a temperature Tmin- The bulk viscosity on the 
small BH background is always larger than the respective one in the large BH background. 
In particular, it can be shown that the T derivative of the quantity C,/ s diverges at Tmin- 
This is the holographic reason for the presence of a peak in C/s near Tc- On the other 
hand, as it is shown in |19], the presence of Tmin (i-e. a small BH branch) is in one-to-one 



correspondence with color confinement at zero T. We arrive thus at the suggestion that in 
a (large N) gauge theory that confines color at zero T, there shall be a rise in C/s near Tc- 
An important ingredient here is the value of the viscosity asymptotically in the small 
BH branch. There the asymptotic value correlates to the IR behavior of the potential. 
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Taking also into account the fact that this asymptotic value is very close to the value of 
the bulk viscosity near Tc, we can derive bounds that suggest that the bulk viscosity cannot 
increase a lot near Tc. 

Drag Force: The drag force calculated from IHQCD has the expected behavior. Al- 
though it increases with temperature, it does so slower than in = 4 SYM, signaling the 
effects of asymptotic freedom. 

Diffusion Time: Based on the drag force calculation the diffusion times can be com- 
puted for a heavy external quark. The numerical values obtained are in agreement with 
phenomenological models |37]. To accommodate for the fact that IHQCD exhibits a phase 
transition around T = 247 MeV (i.e. about 30% higher than in QCD), the results are 



compared using alternative schemes, as proposed in |11]. For example, for an external 
Charm quark of momentum p = 10 GeV (in the alternative scheme) a diffusion time of 
r = 2.6 fm at temperature T = 280 MeV is found. Similarly, for a Bottom quark of 
the same momentum and at the same temperature, r = 6.3 fm. Generally the numbers 



obtained are close to those obtained by [37] and [42|. 



Jet Quenching: The jet quenching parameter of this model, has been also calculated. 



based on the formalism of ||39|, ^ by computing the appropriate light-like Wilson loop. 
q grows with temperature, but slower than the growth of A/" = 4 SYM result. Again 
this can be attributed to the incorporation of asymptotic freedom in IHQCD. Using the 
alternative scheme to compare with experiment, the results are close to the lower quoted 
values of q. For example, for a temperature of T = 290 MeV, which in the alternative 
"energy scheme" corresponds to a temperature of T = 395 MeV in our model, we find that 
q « 2GeV^/fm. 

However, the numbers obtained for this particular definition of jet quenching param- 
eter seem rather low and indicate that this may not be the most appropriate definition 
in the holographic context. There are other ways to define q, in particular using the fluc- 
tuations of the trailing string solution. This is gives a direct and more detailed input in 
the associated Langevin dynamics and captures the asymmetry between longitudinal and 



transverse fluctuations. This can be computed along the lines set in [45, 47, 48 1 and the 
results were recently reported in |]9l|. 
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